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ABSTRACT 


The objective of the research was to develop numerical 
methods of reducing the computational task associated with 
the assignment of damping to the large flexible space struc- 
ture, the dynamics of which, had been approximated to a 
second-order linear matrix ordinary differential equation by 
the finite-element method. 

An efficient decoupling algorithm by which selected 
modes can be damped while the other modes retain their pole 
locations was devised in order to alleviate the computational 
burden caused by the high dimension of the system. Some 
unique properties on a class of the second-order lambda- 
matrices were found and applied to determine a damping matrix 
of the decoupled subsystem in such a way that the damped 
system would have pre-assigned eigenvalues without disturbing 
the stiffness matrix. The resulting system was realized as 
a time -invariant velocity only feedback control system with 
desired poles. Another approach using optimal control theory 
was also applied to the decoupled system in such a way that 
the mode spillover problem could be eliminated. The pro- 
cedures were tested successfully by numerical examples. 

Since the decoupling procedure required only eigenvec- 
tors of the selected modes, the computing time was reduced 
significantly when the number of modes involved in damping 
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assignmeni: was less than one fourth of the total inodes. 
Therefore, in large systems, only a few of the low frequency 
modes need additional damping and the methods of decoupling 
and control developed in this work may be attractive for the 
vibration control of the large flexible space structure. 
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CHAPTER I 


INTRODUCTION 


1 . 1 Overview 

It is well known that one of the problems inherent in 
control of large space structures is that of vibration sup- 
pression. The vibration control is crucial especially when 
the structure is very large and mechanically flexible. In 
particular, space structures such as satellite antennas or 
solar energy collecting panels are expected to be very large. 
Fortunately, with the advent of a space transportation system 
the concept of the large and flexible space structures has 
become more realizable than ever before. However, because of 
the capacity limitations of the space transportation system 
it may be mandatory that the structure be constructed of 
light materials which have very low natural damping. There- 
fore, it is obvious that a certain control action is required 
to provide the large structure with sufficient damping, which 
in turn could maintain the stability of the structure against 
possible disturbances. This comparatively new control pro- 
blem within the field of multivariable control systems was 
considered in this work and a candidate solution to the pro- 
blem was addressed. 

A large space structure may be described as a continue 
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by a set of simultaneous partial differential equations [1]. 
The parameters contained in the equations are, in general, 
continuous functions of the spatial variables. Thus the 
control system is in the realm of a distributed parameter 
system. However, even though the distributed parameter sys- 
tem is rich in theory [2-4] the implementable sensors and 
actuators are difficult to obtain in the infinite dimensional 
space. Consequently, a common approach to modeling is to 
convert the partial differential equations of the distributed 
system into an infinite set of ordinary differential equations 
through spatial discretization [5-7]. A finite number of 
modes are then retained for control. The resulting second- 
order linear matrix differential equation is now of high 
dimension and is difficult to handle on a digital computer. 
Moreover, the dimensionality problem is compounded when the 
second-order matrix differential equation is recast in the 
state variable form. 

The objective of this work is to investigate numerical 
methods based on the properties of the second-order lambda- 
matrices and to devise methods of reducing the computational 
task associated with the finite-element model and the assign- 
ment of damping to the structure. The damping assignment 
will proceed in two different directions ; eigenvalue reloca- 
tion by velocity feedback and optimal control. Each of the 
methods will be based on devising an efficient decoupling 
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technique by which selected inodes can be damped with all other 
modes retaining their pole locations. Although it may not be 
necessary, the theory of the eigenvalue relocation will be 
developed on the basis that tt.e stiffness matrix will not be 
changed. 

In this work we employed the properties of the lambda- 
matrices with the hope that we could alleviate the dimension- 
ality problem mentioned above. Since the inception of the 
term "lambda-matrix" [8] much progress has been made on the 
subject not only in theoretical development but also in the 
application to the system analysis and design. Nevertheless, 
the second-order system which is abundant in the real systems, 
for example, in dynamic structural analysis has often been 
overshadowed by high-order systems. Although the high-order 
system embraces the second-order system as a member, the 
latter may have very unique properties not common to the 
former. This was the primary motivation of this work. 

Therefore, for a class of the second-order lambda- 
matrices, we established some unique theorems which would 
explain the movement of latent roots after a damping matrix 
was added to the undamped system. These theorems were applied 
to compute che damping matrix of the decoupled subsystem in 
such a way that the damped system would have exact pre- 
assigned latent roots without disturbing the stiffness matrix. 
The resulting system was realized as a velocity feedback 


3 


system V7lth desired poles. 

The optimal control theory was also applied to the de- 
coupled subsystem so that the computational load for the 
procedure was reasonable and the mode spillover problem could 
be eliminated. The method is similar to that of reference [9] 
where the modes are decoupled and the optimal control is 
determined for each mode. The main difference is in the 
method of decoupling and the computaional procedure. 

In Chapter II, a large system of second-order ordinary 
differential equations was formulated by the spatial discre- 
tization of the partial differential equations which had been 
assumed to govern the structure. After changing coordinate 
systems the statement of the problem was presented. In 
Chapter III, definitions and theorems on the lambda-matrices 
were collected and followed by the system properties developed 
in this work. Based on this work the decoupling scheme, 
eigenvalue relocation and velocity feedback control synthesis 
schemes were illustrated with examples in Chapter IV. The 
modes decoupling procedure of the state matrix and the optimal 
control of the selected modes were described in Chapter V. 

The conclusion and recommendations for further research were 
given in the final chapter. 

1 . 2 Previous Work and Related Literature 
1*2.1 Control of Large Structures 

In the past several years, diverse groups of scientists 


4 


and engineers devoted much effort on the research of large 
flexible space structures. For a broad view we refer to a 
comprehensive survey paper [II. In addition, we list some 
of the literature which is closely related to our work. 

Most works on the control problem of this d 3 niamic system 
were defined in the state variable form. T.;c primary well- 
developed feedback control techniques with some variations 
wore applied to this problem: they are optimal control [10] 
and pole assignment [11] by constant state feedback. These 
two control schemes were examined in the form of full state 
feedback for possible use in the active vibration control of 
the system and some expected problems were pointed ou: in 
[12,13] . 

Most authors avoided full state feedback for several 
reasons — an awareness of the high dimensionality, the diffi- 
culties associated with the measurement of the states as well 

as the complexity of the Riccati equation. Instead, they 
chose local states or outputs as feedback and/or worked on 
decoupled modes independently. In [14] the local control was 
defined as a control law that included feedback of only those 
state variables that were physically near a particular actu- 
ator. Then a necessary condition for the solution of the 
lineal.* quadratic optimal control problem with the constraint 
of the local state feedback was derived. In [15] the sub- 
optimal output feedback control scheme originated in [16] was 
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applied to the system instead of the computationally difficult 
optimal one [17j. However, this result sacrifices the guar- 
antee of stability for the closed-loop control system. 
Similarly, the inverse optimal control principle [18,19] was 
applied to select the output feedback gain matrix in an 
iterative .manner^ [20] , , In another note [21], a concept of 
member damper controllers was introduced and the problem of 
selecting diagonal velocity feedback gains was formulated 
as an optimal output feedback regulator problem. Finally, an 
optimal control of decoupled modes v/as pursued on the block 
diagonally decomposed subsystems [9,22,23] rather than on the 
full system. Also, the same decoupling technique was applied 
to the control of flexible gyroscopic systems [24,25]. 

For the same reason as in the optimal control, the pole 
assignment was carried out by output feedback to alleviate 
those difficulties mentioned above [26-29] . A direct velocity 
feedback control was suggested as a special case of the output 
feedback [30]. Under some restrictions, even though the 
direct velocity feedback controller can not destabilize any 
part of the system, the exact pole locations can not be pre- 
dicted apriori . Therefore, the initial decision of locating 
and sizing the dampers was based on guesswork and engineering 
judgement. Concerning this problem an idea [31] was pre- 
sented, which enabled the researcher to predict analytically 
the behavior of the closed-loop system by applying root 
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perturbation techniques. However, this was only possible 
when the controller was allowed to modify moderately the 
natural frequencies of the structure. 

Aside from the active control, passive energy dissipation 
mechanisms of the vibrating systems were studied and a future 
design process for sufficient structural damping was sug- 
gested in [32]. Also, a model with the damping matrix pro- 
portional to the positive square root of the stiffness matrix 
was rigorously investigated in [33]. 

1.2.2 Second-Order Dynamic Systems and Lambda -Matrices 

It is well known that undamped linear second-order sys- 
tems possess classical normal modes ,. .whereas in damped 
systems this property is generally violated. However, a 
special class of damped linear systems possesses the classical 
normal modes. A necessary and sufficient condition for this 
class of the systems was shown in [34]. Also, the stability 
was analyzed for second-order systems [35] and for general 
high-order systems [36]. In [37] the computational methods of 
eigenvalues and eigenvectors of second-order systems were 
discussed and these eigen-problems in a user supplied inter- 
val were examined in [38] . Finally, controllability and 
observability [39], a solution of the eigenvalue problem [40], 
and a modal analysis for the response of this system [41] were 
studied on linear gyroscopic systems. 

On the other hand, the ana'vaig of vibrating systems by 
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lambda-matrices was introduced in earlier works [8,42]. 

This mathematical area was rigorously developed by a group 
of mathematicians and, as a result, a comprehensive treat- 
ment was provided in [43]. Along with this development, the 
algebraic theory of the matrix polynomials [44] , algorithms 
of their solvents [45,46] and spectral factors [46,47] of the 
corresponding lambda -matrices were established. Also, trans- 
formations of solvents and spectral factors were developed 
in [48], Consequently, as multivariable theories attracted 
much attention from engineers lately, those lambda -matrix 
related theories began to be applied to such engineering 
fields as filter design [49,50], partial fraction expansions 
[51] and multivariable control systems in general [52,53]. 
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CHAPTER II 


FORMULATION OF THE PROBLEM 

2. 1 Description of the Mathematical Model 

The large space structure is assumed to be described by 
the following system of partial differential equations [1,5]; 


M (p)3iu(E^ 


at' 


^ at 


+ KQU(p,t) =fQ(p,t), (2.1) 


where u(p,t) is the displacement of an arbitrary point p of 
the domain D off its equiribrium position, Mq(p) is the 
distributed mass, Kq is the time -invariant symmetric non- 
negative differential operator of order 2P, and fQ(p,t) is 
the distributed control vector. The damping term Cg— 
which is S 3 rmmetric and represents the internal structural 
damping, is thought to provide the structure with very weak 
mode damping. The displacement u(p,t) must satisfy the 
boundary conditions at every point of the boundary S of the 
domain D; 

u(p,t) =0, i = 1,...,P (2.2) 

where B^, i=l,...,P are linear differential operators. 

Because of the theoretical difficulties and practical 
complexities arising from the implementation of the infinite 
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dimensional control system [54] we approximate the infinite 
dimensional system by a spatial discretization of the partial 
differential equations. Among a variety of available methods 
we choose a finite- element method such as Galerkin's method 
[55], 

It is assumed in the method that an approximate solution 

A 

of (2.1), u(p,t), which also satisfies the associated boundary 
conditions (2.2) could be expressed as 

u(p,t) = I <{..(p)v.(t) (2.3) 

i=l ^ ^ 

where (f)^(p) , i=l,...,n are comparision functions^ depending 
only on the spatial coordinates and v^(t) , i=l,,,.,n are 
time “dependent generalized coordinates. Since (2,3) is only 
an approximate solution it may not satisfy (2,1) exactly. 
Hence, in Galerkin's approach, the v^(t) is chosen to mini- 
mize the mean-square equation error when (2.3) is substituted 
into (2.1). 

Once v^(t) is determined, substituting (2.3) into (2,1), 
pre-mult ip lying both sides by 'l’^(p) and integrating over 

^ Comparision functions are any arbitrary functions satisfying 
all the boundary conditions of the eigenvalue problem and are 
2P times differentiable over domain D[6, P, 140]. In 
Galerkin's method they could be linear combinations of piece- 
wise linear functions or cubic splines, etc. 
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OF POUR QUALITY 

the domain D lead to 

(j)^(p)MQ(p)<j)^(p)dD + I v^(t)[ <|)^(p)C.tj).(p)dD 

i— 1 D i=l ■' D ^ u 1 . 

+ I (J>J(p)kQ(i) ,(p)dD = [ <()^(p)f.(p,t)dD,.=l n 

1=1 JD 

(2.4) 

where the superscript t denotes the transpose of the matrix. 
Thus , by introducing the notations 

f ^Cp)Mj(p).|.j^(p)dD - 

j^*J(P>Co*i(p)dD - C^.. 

■^♦J(p)KQ*3^(p)dD - K^. . 

. ocj)^(p)fo(p,t)dD = f^(t) , 

(2.4) reduces to 




( 2 . 5 ) 


/N 

where fj,(t) is the generalized forces associated with the 
generalized displacement v^(t) ; It must be noted that 

because of the charateristics of Mq(p) 


i 

1 

t: 
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and Kq. Now, by letting 

v(t) = [Vj^(t) , . . . ,v^(t) 

and 


f(t) = [f^(t) , . . . ,f^(t)]^, 

(2.5) can be written in the matrix from 

Mv(t) + Cv(t) + Kv(t) = f(t), (2.6) 

where the mass matrix is S3nnmetric positive definite, 

the damping matrix is symmetric, and the stiffness 

matrix KeR’^^^ is symmetric semi-positive definite. 

By changing coordinates v(t) to x(t) by x(t)= M'^v(t) 
and pre-multiplying both sides of (2.6) by M we have 

Ix(t) + Cx(t) + Kx(t) = f(t) (2.7) 

where I is the (nxn) identity matrix, C = K = 

and f(t) = M ^f(t). The change of the coordinates preserves 
the system eigenvalues as well as the symmetry of the system. 
The associated homogeneous system of (2.7), 


Ix(t) + Cx(t) + Kx(t) = 0 


( 2 . 8 ) 


can be transformed to a corresponding lambda-matrix by 
assuming a solution x(t) = Xq exp (At) or by taking the 
Laplace transform of (2.8) with zero initial conditions. 
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The lambda-matrix thus obtained is expressed as 


A(X) = + CX + K, (2.9) 

and this monic S3nnmetric second-order lambda-matrix will be 
the main subject of chapter III. 

2.2 Statement of the Problem 

Consider for a monent a single degree of freedom 
spring-mass -damper system with free vibration. This system 
is characterized by a homogeneous differential equation 

X(t) + 2?to^x(t) + = 0, (2.10) 

where x(b) denotes the displacement from the equilibrium 
position, ? is known as viscous damping factor and is 
the natural frequency of the system. Substituting the 
solution x(t) =aexp(Xt) into (2.10) yields the character- 
istic equation, 

X^ + 2^ojj^X + 0)^ = 0, 

with two roots, X, = -Cw + jiu /l-?^ and X« = - jw /l-?^ 

where j=/-l. Hence, the solution of (2.10) can be written as 

X(t) = aj^exp(Xjt) + a2exp(X2t), (2.11) 

where a, and a„ are constants depending on the initial dis- 



placement and velocity. 

Now, from the above elementary solution of the homo- 
geneous system (2.10) we can observe several characteristics 
of the system. First of all, the behavior of the system 
with non-zero initial conditions (disturbance) depends solely 
on the value of the viscous damping factor; When ^>1 the 
motion of x decreases monotonically with increasing time. 

This is the overdamped case. When 0<?<1 the system is under- 
damped and it decays sinusoidally. In case ? = 1 the system 
is critically damped and is the limit between the regions. 
When ?<0 the system becomes unstable. So, it explodes(5<0) 
or oscillates (t=0) under the presence of a small disturbance. 
Therefore, the stability of this single degree of freedom 
system is determined only by the sign of it is stable if 
?>0, otherwise it is unstable. 

Furthermore, the damping factor appears explicitly in 
the equation so that the designer can choose this quantity 
before testing and can predict the system response. This 
property is especially desirable when the design specifica- 
tions are very stringent as in the case of the large space 
structure. Also, this arbitrary assignability of damping 
may be critical under the system parameter uncertainty, in 
which case the margin of stability plays a great role in the 
design of robust system. 

Not surprisingly, none of these properties of the scalar 


system belongs to the coupled multi-degree of freedom system 
as in the example of the coupled two-degree of freedom system 
shown below[35] ; 


1 0 


8 0 


24.8 14.5647 


x(t) + 


x(t) + 


_0 1 _ 


_0 -2_ 


J4.5647 9.2 _ 


( 2 . 12 ) 

The eigenvalues of this homogeneous system are Aj 2 ” 
and ^ = -2±2j, so that the system is asymptotically 
stable. Nevertheless, the damping matrix is not positive 
definite, which indicates a completely different property 
from that of the single degree of freedom system. 

Concerned with the development of unique properties 
on a class of second-order djmamic systems and with damping 
assignment to the systems we define the problem as follows: 
i) to develop properties of a homogeneous second-order 
symmetric system 

Ix(t) + Cx(t) + Kx(t) = 0 (2.13) 

where is symmetric and is symmetric semi-positive 

definite, ii) to find a s 3 mimetric damping matrix CeR^^^ with 
given K such that (2.13) has some (or all) pre-assigned system 

eigenvalues, iii) to determine a collocated velocity only feed- 
back control vector f(t)eR™^^ of Ix(t) + Kx(t) = BQf (t) , 
so that the closed-loop system is identical to the system 
found in ii) , where the output y(t)eR^^ is defined by 
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y(t) » BQX(t) , 


and control vector f(t)eR™^^ by f(t> = Gy(t) with a time- 
invariant gain matrix GeR™^. 

In the second part of the work, we define the system 
equation in state space form: 


n(t) = An(t) + Bf(t), BsR^^^°* 


(2.14) 


where'’ 


n(t) = 


x(t) 

x(t) 





We seek a control vector f(t) which is a function of state 
vector? f(t) = rn(t), TeR’^^^ and minimizes a certain quad- 
ratic cost functional with constraints (2.14), but does not 
change all of the eigenvalues of the closed-loop . system 
matrix except for a few modes selected for control. 


^ The notations 0 and I (or 0 and I ^) indicate nxn 

n n nxm nxm 

zero matrix and identity matrix (or nxm) . These will be 
used throughout the work unless otherwise stated. 
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CHAPTER III 


LAMBDA-MATRICES AND LINEAR SYSTEMS 


3 . 1 Lambda-Matrices, Matrix Polynomials and Solveata 

It was shown in the previous chapter that a second-order 
lambda-r trix, 


A(X) 


+ CX + K 


(3.1) 


was obtained from the Laplace transformation of the system 
equations (2.8) with zero initial conditions. In this 
section definitions and theorems on the lambda-matrices, 
matrix pol 3 niomials and solvents, which are essential to the 
development and design of damping system, will b^ summarized 
for the continuity of the presentation. Proofs of the 
theorems and more rigorous treatment of the general case 
are given in Appendix A for lambda-matrices and in Appendix B 
for their solvents. 

A latent root of a lambda -matrix A(X) is defined as a 

nxl 

number XeC such that det A(X)=0 and a non-zero vector yeC 
is called a right latent vector corresponding to the latent 
root X if A(X)y=0^^2 where det denotes the determinant. 
Similarly zeC’^^^ is a left latent vector if ( X ) z=0^2j,^ • 

As in the state space modeling of any matrix differen- 
tial equation the lambda-matrix / (X) can be associated with 
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the block companion matrix A eR which has the following 

structure; 



It is shown in Appendix A that the 2n ei.genvalues of Aq are 
the latent roots of A(A) counting multiplicity. Therefore, 

the terms "system eigenvalues" and "latent roots" of A(A) 
will be used interchangeably throughout this work. 

Associated with the lambda-matrix (3.1) are two types 
of polynomials, a right matrix pol 3 niomial, 

A^(X) = + CX + K 

and a left matrix pol}niomial, 

Al<X) == X^ + XC + K 

where XeC’^^^. If "^^^^“®nxn 

arr called a right and a left solvent of A(X) respectively. 
The right (left) solvent is viewed as a matrix root of the 
right(left) matrix polynomial and the solvents retain com- 
plete information of latent roots and latent vectors. These 
characteristics of the matrix pol 3 momials are explained- in 
the following theorems . 

Theorem 3.1 . If A(X) has n linearly independent right 
latent vectors y^, i=l,...,n corresponding to the latent 
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roots of A(X), (X^ is not necessarily distinct) 

-1 

then YAY is a right solvent, where Y == [y 2 *«**>y^] 

A = diag(X;i^, . . . , X^) . 

Theorem 3.2 . If A(X) has n linearly independent left 
latent vectors Zj^, i^lj.-.jn corresponding to the latent 
roots of A(X) , X^, i=l,...,n (X^ is not necessarily dis- 
tinct) then Z~^kZ is a left solvent, w>.ere Z = [z^,...,z^] 
and A = diagCX^ > . • . »X^) . 

Proofs of Theorems 3.1 and 3.2 are included in Appendix 
B. Notice that these theorems are special cases of Corol- 
laries B.5 and B.6 in the appendix, respectively. 

Theorem 3.3 . A(X) is divisible on the right(left) by 
IX-R (IX-L) if and only if Pv(L) is a right (left) solvent 
of A(X) . 

Proof; See Appendix C. 

3.2. Properties of A(X) = IX^ + CX + K 

Definitions and theorems in the previous section can be 
directly extended to n-th order lambda -matrices , whereas the 
properties developed in this section may be valid only on 
the second-order real symmetric system. However, these 
properties will serve as the foundations for the damping 
system synthesis which will be developed in the next chapter. 

Theorem 3.4 . Suppose all latent roots of A(X) exist in 
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distinct complex conjugate pairs. If a constant matrix 
R^C^xn ^ right solvent of A(X) and R has one latent root 
and corresponding latent vector out of each complex conjugate 
pair , then 

i) R* is a left solvent of A(A) 

ii) A(A) = (IX - R*) (IX - R) (3.2) 

iii) R*R = K (3.3) 

iv) R* + R = -C (3.4) 

where * indicates complex conjugate transpose. 

Proof; i) According to Theorem 3.3 A(x) can be factored 
on the right by (IX - R) , i.e., 

A(X) = Q(X) (IX - R), 

where Q(X) is a quotient lambda-matrix. Upon taking the 
complex conjugate transpose of A(X) it follows that 

A*(x) = (IX* - R*) Q*(X). 

But, since C and K are real symmetric matrices, A*(X)=A(X*) . 
Therefore, 

A(X)'= [A*(X)1* = A*(X*) = (IX - R*)Q*(X*) 

Thus, R* is a left solvent of A(x) by Theorem 3.3. 
ii) The factorization is obvious from the fact that R and 
R* have a disjoint and exhaustive set of latent roots and 
Q(X) is of first order 
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iii and iv) Equating the coefficient matrices of (3.1) and 
(3.2) completes the proof. VVV 

Lemma 3 .5 [58] . Suppose there exists a matrix U-stich that- 
A -h 

U = RK . Then, under the same assumption of Theorem 3.4 


i) U is unitary 

ii) n iX, ( " = n u? (3.5) 

i“l' i=l 1 


where and w? are eigenvalues of R and K, respectively 

and lx. I denotes the absolute value of X . . 

' l' X 

Proof; i) UU* = RK'^R* 

= RR“^(R“^)*R* {from (3.3)} 

= I 

U*IJ = k"^r*rk’^ 

- K“^KK~^ {from (3.3)} 

= I 

ii) From i) we have R=UK^. Taking the determinant and 
absolute value of this equation in order, we have 


or 


I det R ( = 1 det U | • | det ( 



n 

H ujj 

i==l 


where positive definiteness of K was utilized. VVV 

•This lemma establishes the relationship between natural 
frequencies and latent roots of the damped system and will 
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play an important role in assigi'.ing arbitrary damping, which 
will be shoxTn in the next chapter. 

Theorem 3.6 . Let y^, i=l,...,n be latent vectors of A(X) 
with y£*y£ “ 1 and define 3^ = y^ Ky^^^. Then, 


(0^ . < 3 • < 

min “ - 


w. 


max’ 


i=l , . . . ,n^ 


(3.6) 


where . and are the minimum and the maximum eigen- 

min max ® 

values of K, respectively. 

x^Kx 

Proof: Let the Rayleigh's quotient f(x) = x ■ where 

X ^ I 

xeC . When x is replaced by eigenvectors of K the 
Rayleigh's principle provides the proof immediately. 

Tlxl 

Lemma 3.7 . Let X^eC and y^sC , i=l,...,n be latent 
roots and a set of latent vectors of A(X) such that y^*y^=l , 
and define a. = %y.*Cy. and 3. = y.*Ky., i=l,...,n. Then, 
the system can be classified as follows: 

i) When ct^>0 and a?>3j|^ 

a) X^ = -a. ± /a? - 3^, 

b) xtxT = 3., where xt = -a. + /a? - 3., 

= -a. - /oc^. - 6^ and 

c) the i-th mode is overdamped. 

ii) When a^>0 and 

a) X^ = -a^ ± j/3^ - 
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b) and 

c) the i“th mode is underdamped. 

iii) When = 0 the i-th mode is oscillating, 

iv) When < 0 the i-th mode is unstable. 

Proof : Since and are latent roots and latent vectors 

of A(X) respectively, + K)y^ == 0. After .pre- 

multiplying y^* xc follows that 

X? + 2a.A. + 3- “ 0, i=l,...,n. 

Therefore, the latent roots can be expressed in terms of 

and 3^, i.e., = -a^ ± /a| - 3^» which provides the proof. 

VVV 

Lemma 3.7 reveals that the positive definiteness of 
the damping matrix C is not a necessary condition for the 
stability of the system. In other words, even though the 
matrix C is indefinite it is still possible that there are 
some latent vectors y^ such that y^*Cy^ = a^>0, i='l,...,n, 
which is a necessary and sufficient condition for a stable 
system. The result obtained in (2.12) can now be explained 
by this theorem. 

Also, Lemma 3.7 combined with Theorem 3.6 shows the 
feasible region of latent roots of a damped system in terms 
of natural frequencies of the undamped system. After a 
damping matrix G is added to the undamped system, 

A'(A) = lA^ + K, the new latent roots are bounded by the 
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concentric circles with radii of maximtim and minimum natural 
frequencies of the undamped system. This is illustrated in 
Fig. 1 and the system A(X) = + Cx + K is classified 

according to the asymptotic stability of the modes in Fig. 2. 
Furthermore, Lemma 3,5 indicates that the new latent roots 
are totally governed by the natural frequencies of the 
undamped system. However, their precise locations are 
affected by the new latent vectors which will be changed 
after damping is introduced. The shift of latent roots 
due to damping will be explained after Lemma 3.8 (classical 
normal mode) is presented . 

Lemma 3.8 . If the damping matrix C has the same modal 
matrix as the stiffness matrix K 

i) jX. I ^ = 01 ? , i=l, . . . ,n when oi? > a? and 
ii) ~ 0 )^ , i=l,.,.,n when w? < , 

X 1 1 ‘ XI 

where oi? and 2a^ , . is=l, . . . , n are eigenvalues of K and C, 

+ - 

respectively, and X^ and X^ are the same as in Theorem 3.7. 

Proof: By the hypothesis, let y,eC^^^ in Lemma 3.7 be a 

normalized eigenvector of both matrices C and K. Then, 

Bi = w?, i=l,...,n. Substitution of w? for 8^ in Lemma 
3.7 gives the result of Lemma 3.8 directly. VVV 

From this lemma we see the movement of latent roots 
after damping is reinforced in case of the classical normal 
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ZiSH^eJ,: Feasible region of latent roots of the 

underdamped system A(X) = 


+ CX + K. 




X 

Stable & 
over damped 


Oscillatory 


- 03 . 


max 


2 

Figure 2 ; Classification of system A(X) = lA +CA+K 
by the stability of the modes. 












modes. In this case the absolute value or the modulus of 
is equal to the eigenvalue of K, whereas in the general case 
as seen in Lemma 3,5 multiplication of the moduli of x^ is 
equal to the multiplication of the eigenvalues of K. The 
shifts of latent roots for these two cases have been illus- 
trated in Figs . 3 and 4 . 

3.3 Bounds on Latent Roots of A(X) = IX^ + Cx K 

Careful examinations of the solvents of the linear 
second-order system unveil some interesting relations between 
solvent and coefficient matrices of the system. Based on 
these relationships some bounds on the eigenvalues of the 
damped system will be developed in this section. We separate 
solvent into the real part and the imaginary 

part R^eR^^^. Now, from Theorem 3.4 it can be shown that 
Rj is symmetric and the following relations hold for the 
system; 


4 -^ 

(3.7) 

R^ Rj. = RjRj^ 

(3.8) 

Rj^ + Rp_ = -G. 

(3.9) 


Before developing the bounds on the eigenvalues of the 
system we introduce theorems necessary for the development. 

Theorem 3.9 . Let f(x) * ^ for xeR^^^, x 0 and 

XX 
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_s_££_: Shifts of latent roots after damping 
assigned(Classioal normal mode case - C and K h 
Che same eigenvector matrix). 






If xeR^^^ is an eigenvector of %(A + A^) , then f(x) 
is the eigenvalue of %(A + A^) corresponding to x. 

A A 

Proof: Let \ and x be an eigenvalue and corresponding 

eigenvector of %(A + A ). Then, 


r _ x'^CA + A^)x _ 1 / x^Ax , x^A^x ^ x^'Ax _ ^/''x 

A = z — = t( — p — T - — r — ) = — z — “ f(x) 

-AUA Z aCa aCa aCa ' ' 


— 

2x X 


Corollary 3.10 . Let f(x) 


XX XX 

.t^ 


X X 


VVV. 


X 


Rj^j 


,nxl 


for xeR X 0. If 


x^x 

^ nxl ^ 

xeR is an eigenvector of (-C) then f(x) is the eigenvalue 

A 

of (-%C) corresponding to x, where C is the damping matrix 
of A(A). 


Theorem 3.11 [59] . Decomposing an arbitrary matrix A into 
A = + jH 2 » where and H 2 are Hermitian^ ; then for every 

eigenvalue A of A we have 

i - ^max<H2) 


where A A„„„, Re, and Im denote the minimum, the maximum 
min max 

eigenvalue, real part, and imaginary part respectively. 
Proof; See Appendix D V7V 

^ Note that and H 2 are defined as = %(A + A*) , 

^2 ^ - A*) and that in general they are not real. 
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Lemma 3.12 . Let be a right solvent of A(X) . Then 

the real parts of the eigenvalues of R have the fallowing 
bounds ; 




< Re X(R) < 


where C is the damping matrix of A(A). 

Proof; Since Hi = %(R+R*) “ -%C from (3.4) Theorem 3.11 
provides the proof immediately. VVV 

This lemma provides the bounds on the real parts of the 
eigenvalues of the damped system. 



CHAPTER IV 


SYNTHESIS OF DAMPED SYSTEM BY EIGENVALUE RELOCATION 

4.1 Decoupling of the Stiffness Matrix 

One of the major problems associated with the large 
scale system is that of high dimensionality. One possible 
approach to reducing the dimensionality is to decouple the 
system equations into subsystem equations where the col- 
lection of subsystems retains the eigenvalues of the original 
system. The most obvious method of decoupling the large 
system is to compute the eigenvectors of the system and then 
reduce the system equations to those of the Jordan form, 
a total decoupling. It is not necessary and probably not 
preferrable to decouple completely but to decouple only those 
modes that are to have additional damping. This partial de- 
coupling has an advantage in that the feasible locations of 
new eigenvalues are more flexible than the classical normal 
mode case as shown in Figs. 3 and 4 of Chapter III. In this 
section the theoretical background for the partial decoupling 
will be developed. 

Theorem 4.1 . Let D^eC’^i^’^i, i=l,...,A consist of Jordan 
blocks of KeR^^’^ in such a way that no Jordan block is split 
among D^ s and ^£^n^» n. If i~l i exist for arbit- 
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rarily chosen matrix K can be -block- diagonalized 

through a similarity transformation to of the form, 

“ ®diag^%l’ • * • 

where indicates a block diagonal matrix and 

“ Wi’'- i 1 

Proof: Let ^ be an eigenvector matrix of K, 'I'g = 

(f 2 _» • • • and ’ * * ' * Sirice exists from 

the existence of i=l,..,,5, Kg can be expressed as Kg= 

equality D = $“^K$ is substituted, Kg= 

TKT ^ is obtained where T = Since T~^ exists, 

the proof of the theorem is completed. VVV 


The matrix in the above theorem is arbitrary and 
T can be computed directly. However, the decoupled block 
Kg., is, in general, not symmetric. The following lemma shows 
a way to obtain symmetric blocks. 


Lemma 4.2 . Let KeR^^ be symmetric and of simple structure 


and D^eR^i^^i be a diagonal matrix which consists of eigenval- 

a 

ues of i-th block with .S_n.=n. Also let $ be an orthogonal 

IL 1. X 

eigenvector matrix of K. Then, K can be block diagonalized 


into KgeR' 


nxn 


by a similarity transformation T such as 
Kg = TKT , where Kg = ^diag ^^11^1^11 ’ 


33 


“Si 

OF POOR QUALITY 


T = 




♦ll*U 

« 

« 


I is appropriately block partitioned as 


? = 




# 


1 ^jll ••• 

and i=l,...,H are orthogonalized matrices of ^..eR’^i^^i. 

Proof; Substituting for i=l,...,£ in Theorem 4.1 
gives Kg^ ^ii^i^ii 


T 



^11^11 • • • ^ll^U 

• • 

• • 

• • 


vvv 


If the interest is in assigning damping to a few of the 
low frequency modes, it is not necessary to compute all of 
the eigenvectors of K. Since most of the energy of the 
system is in the lower frequency modes, a reasonably con- 
trolled structure would require control of the first few 
modes with the exception of the rigid modes. Under these 
circumstances the following ^•’^eorem is useful for the eigen- 
value assignment; 


diag^^q ’ 


-In-q) . S = km 


-1 


and 
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Theorem 4.3. 


Let E = B 


OF POOR QUA JTY 


T = S + %E where is an eigenvector matrix of KeR^^^. 

Then, T will block diagonalize K into the form Kg = 

-1 

(Kgi , Kg2^ under a similarity transformation; Kg = TKT , 

where Kg^sR*^^*^ has eigenvalues corresponding to the first q 

eigenvectors of and the rest of them belong to Kg 2 e 
P^(n-q)x(n-q) ^ 


Proof [60]: Let $ be partitioned into four blocks 

“^ll ■‘^12 
^ ^21 ‘^22 

where ^ud q)x(n q) ^ Then, with E as defined, 

T = S + 

= 2[$E + 

‘*’11 ^qx(n-q) , 

= 0 

°(n-q)xq "***22 

Since K = $D$ 

Kg = TKT'^ 

" ®diag^‘*’lA‘*’ll' '*’22°2‘*’22^ 

provided that no eigenvalue is common to both and D 2 . 
when we denote K^. = t() . ,D. cj) . K i=l,2, the proof is completed. 

vvv 


35 


ORfGf.Vv? : 
POOR Q 




V f 
■J 


'iY 


The similarity transformation in Theorem 4 . 3 can be 
obtained from the knowledge of the eigenvectors corresponding 
to those modes that are to be decoupled. Let be the 
eigenprojector^ of the q modes that are to be decoupled; 



$ 



^(n-l)xq 


0 

0 


qx(n-q) 

n-q 


<5 


t 


♦ll*u ♦u'*'21 
♦21*11 *21*21 


-1 t 

where it is assumed that is normalized such that $ =?$ . 

The matrix S defined earlier is then 



and thus 


T = Pj3 + 


qx(n-q) 


0/ X 1 

(n-q)xq n-q 


Since the matrix contains only the first q eigenvectors 
t t t 

of K, ^21^ ' been shown to be constructed from 

those q eigenvectors. 

^ The eigenproj ector corresponding to the first q eigen- 

^ t 

values of a simple matrix K is defined as Pjj = 

where (j)^ is an i-th normalized eigenvector of K. For the 

general case and properties of the eigenproj ectors, see [51]. 
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The computational procedure given above can be utilized 
to find which contains the desired eigenvalues that are 
to be moved av^ay from the jw axis of the undamped system. 

The computational technique to assign the eigenvalues will 
be discussed in the next section. 


4 . 2 Computational Procedure for the Damping Matrix Determination 
Assume that the system stiffness matrix K is block 
diagonalized into and that tias the eigen- 

values of the system that are to be moved. The subsystem 
matrix can now be further decoupled and it wxll also be 
assumed that is block diagonalized to (2x2) matrices 
Ktj . with two distinct real eigenvalues^ . The theory given 
previously can now be utilized to assign the eigenvalues 
of the low-order subsystems. 

Let and y^eC for i « 1,2 be the latent roots 

and a set of independent latent vectors of the desired sub- 

* 2x1 

system with y^y^^ = 1. Define a new set of vectors w^eC 

such that Wj^ = 0^.y. for i = 1,2 where <{).. is a normalized 

eigenvector matrix of Kgj . Let each component of w^ be 

expressed in polar form. 


’^ik = ^ik 


exp(j0>^j^) , i,k = 1,2 


(4.1) 


^ Subsystem matrix Kgj has real eigenvalues when the stiff- 
ness matrix K is real symmetric and of simple structure. 
This can be justified by Lemma 4.2. . , 
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where and 0^^^ are the modulus and argument of the complex 
number w^j^, respectively. Since 7^*7^ ~ 1 then w^*w^ = 1 and 


4i + 


2 2 
^21 ^22 


1 

1 . 


Furtheinnore, 


(4.2a) 

(4.2b) 


”!2”U 

+ 

W*jWj2 

= 2auai2=°=01 

(4.3a) 

^22^21 

+ 

”2l”22 

“ ^®21®22®°“®2 

(4.3b) 


where 0^ - B ^2 " ®2 ° ®21 " ®22’ 

2x2 

On the other hand, assuming. R. eC is a solvent of the 

* . *3 

damped subsystem with stiffness matrix Kg ^ , from Theorem 3.4 
we have 

= R^R. . (4.4) 

3 3 

2 2 

Now let W = [wj^ W 2 ] , Y = [y^ y 2 ] , = diag(u^, ( 02 ) and 

2 

A = diag(Xj^, X 2 ) where and X^, i=l,2 are eigenvalues of 

K-n . and latent roots of the damped subsystem. Then from the 
«3 

relations, Kg. = K- = YAY"^ and Y = we have 

333333 33 

W*A„W = A*W*WA . 

K. 

This matrix equation is then written. as follows; 


LHS = 


2 2 
“ 1^11 


X 2^2 
+ tOoS 


2“12 


“Hit'll + “2''22"l2 


“l"tl”21 + “2'”!2”22 


2 2 2 2 
‘^ 1^21 ^" 2^22 
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RHS 




^ 1^2 ^^ 11^21 ^ 12 ^ 22 ^ 


^2^1 ^''21^11 ^ 22 ^ 12 ^ 


^2h 


Equating elements of (4.5) and noting that the matrix is 
symmetric we have three equations : 




(4.6a) 


2 2.2 2 I , 

^ 1^21 * ^ 2^22 ’^ 2 ^ 


(4*6b) 


wflW2l(X{\2 - “ Wf2”22'“2 


XfXj) 


(4 .6c) 


Since '^*2^22 ^ ^ complex niimber s can be defined from 

(4.6c) as 


4 

s = 


^1 


“2 - Hh 


w 


!2”: 


22 


^1^2 “ “1 


(4.7) 


The number s can also be obtained from (4.1) as 


^ By the independency of w^ and W 2 both w ^2 ^22 

Wii and W 2 j^)- cannot be zeroes at the same time. If either 
w ^2 ^22 is zero then ^^ 2 ^ 2 ^ = 0 except in two cases, i.e 

W = [q 2 ^ and W = [2 • It now follows from W = 

that Y = [02 ‘?>2^ ^ f‘J’2 ’ respectively. These are 

the classical normal mode cases. 




- ^ 


QR\ami. PPM' 

OF POOR QUAUT/ 


^12^22 


(cos0 + jainO) 


where 0 = 0^ + ©g “ tan"^(Im. s/Re s) 


(4.8a) 

(4.8b) 


Theorem 3.4 also gives 


Sj =■ R? + 


(4.9) 


and by a similar substitution as in (4.5) 


" -(W*)‘^A*W* - WAW"^ (4.10) 

A = 

or with C„ . = (i..C„.d).. 

Bj ^3 33 

A*W*W + W*WA = -W*CgjW. (4.11) 


Substituting the previously defined variables into (4.11) 
we have a set of nonlinear equations; 


2 

ail 

^^11^12^°®®! ^12 


~ — 1 

A 

n 


2(^2 

2 

^21 

2 

^^21^22^^°®®2 ^22 


®11 

A 


202 

Re s 

^21 ^11 
~^cos0o4':~ 1 

^22 2 ^12 ^ 


=12 


Re t 

Im s 

^21 . ^11 
-^sinO^-K^sinOT 0 

^22 ^12 ^ 


C 22 


Im t 


where , i,j=l,2 are elements of Cgj , i=l,2 

parts of and t = (A^ + A 2 ) (s + 1) . 


(4.12) 


are real 
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The computational procedures for the subsystem can be 
summarized as follows: 

2 

1) Compute the eigenvalues and normalized eigen- 
vectors i = 1,2 of Kg j . 

^2 2 

2) Select new eigenvalues and compute = |X^| , 
i = 1,2. 

3) Solve the set of equations for a|j , i,j = 1,2: 



4) Compute Re s and Im s from (4.7) 

5) Compute 0 from (4.8b) 

6) Solve the set of nonlinear equations given in (4.12) 

A 

for the coefficents C.., i,j = 1,2. 

7) Compute the damping matrix C„ , = 

J 33 333 


The subsystem is then transformed back with the other 
subsystems to obtain the damped system equation. 


4 . 3 Velocity Feedback Control Scheme 

In the two previous sections it was shown that the 
stiffness matrix K can be decoupled into block diagonal 
matrices , K 


Tt 


= (j) . ,D.0 . . , j = through orthogonal 

^ J J J J J «J 

similarity transformation T = 'i'g?^ (see Lemma 4.2) and t 
sufficient damping can be assigned to each subsystem by 
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shifting system eigenvalues. In this section, damping as- 
signment will be realized through a velocity only feedback 
control. For simplicity of the implementation, collocated 
control (in which the actuators and sensors are placed 
together) will be employed. 

Neglecting the low natural damping the large flexible 
system of which the vibration is to be controlled is express- 
ed as 

X + Kx = Bgf 

y = B^x (4.13) 


where xeR^^^, yeR^^^ and BqsR’^^™ are generalized 

displacement, control input, output, and actuator matrix, 
respectively^ . Notice that the form of sensor matrix Bg of 
the collocated control. With the state vector n defined by 

r *1 


n 


A 


X 


(4.14) 


the system described in (4.13) is now written in the state 
space form as 



6 

I 


~ — 
0 

^ = 

n 

n 

n + 

nxm 

-K 

0 

n_ 


Jo _ 

y = 

mxn 


n . 



(4.15) 


^ The notation of vector y is different from that in Chapters 
II and III. This notation will be used in this chapter only. 
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Let be the time- invariant velocity feedback gain 

matrix, then this control law will be 

* = bJ] n 


Therefore, the closed loop system is described as 


n 






(4.16) 


and its eigenvalues are identical to the latent roots of the 
corresponding lambda -matrix , 

A(X) « IX^ + CX + K 

where C = -BqGBq. (4.17) 

Thus the remaining problem is to decompose C into 
BqGBq, so that resulting feedback gain matrix G together 
with a sensor matrix Bq and an actuator matrix Bq will pro- 
vide the system with sufficient damping by means of shifting 
system eigenvalues to pre-assigned locations. 

We will now show this prodecure in a constructive way. 
Suppose that -the stiffness matrix is decoupled into three 
subsystems Kg 2 

that Kg^ and Kg 2 are (2x2) subsystems whose modes are to be 
damped and Kgg has the rest of the -modes which will remain 
unchanged. Then, according to Lemma 4.2 


Kgg by Lemma 4.2. It is also assumed 
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where 


Kg = TKT 


-1 


OF POOR QUALITY 


and 




T = 




0 

0 


hi ° 


^22hh2 0 


^ 33 ^ 3*^33 


0 


0 <p22 

0 0 ^23 


-t 

$ , 


(4.18) 


On the other hand, for subsystems Kg ^ , j =1,2 damping 
matrices Cgj , j =1,2 could be determined by the procedures 
described in Section 4.2, and they are assvimed to be written 
as 



• - Cn . d) « • 


j = 1, 2 


(4.19) 


Therefore, in the generalized coordinate system the damping 
matrix will be 
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Substituting (4.18) and (4.19) into (4.20) we have 


C 


hi 

^12 

hi 

hz 

hi 

^32 


'hi » 


^11 ^21 ^31 

^ S2 

L_ _ 


$^2 ^^2 ^^2 
— . 


(4.21) 


Thus, by comparing (4.17) with (4.21), a solution to the 
collocated velocity feedback control problem is as follows; 


Gain matrix; G = ^diag^ ”^B1’ 
Actuator matrix; Bg « c {)2 
Sensor matrix; Bg 


■^B2'" 


•• ‘*’2k^ 



) 


where k is the number of (2x2) subsystems to which damping 
is to be assigned and is a normalized eigenvector corre- 
sponding to the i-th mode. 


4.4 Illustrative Example 
Example 1 ; 


The example chosen for computational purposes has a 


mass matrix I and a stiffness matrix K 


9 

-5 


K = 


-5 

11 -6 
-6 13 -7 

-7 15 -8 

-8 17 

-9 




-9 

19 
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which has natural frequencies given in Table 1. The zeroes 
of the matrix K have been deleted for the sake of brevity. 
The two lowest modes are to be damped with the real parts 
of the damped eigenvalues at -0.5 and -1.0. 


Table 1. Natural Frequencies 

in Example 1 

Mode i 

(Oi 

1 

1.135 

2 

2.215 

3 

3.175 

4 

3.980 

5 

4.687 

6 

5.470 


Through the computational procedures given in the pre- 
vious section the velocity only feedback gain matrix G and 
the force actuator matrix Bq (the sensor matrix Bq assuming 
the sensors and the actuators are collocated) were obtained 
as follows: 


G = 


1.053 

-1.099 


-1.099 

1.949 
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The damping matrix C and the eigenvalues of the closed-loop 
system were computed and listed in Table 2 and Table 3, 
respectively. The eigenvalues of the closed-loop system show 
that the two lowest modes get the pre-as signed amount of 
damping exactly while the rest of the modes remain unchanged. 


Table 2. The Computed Damping Matrix C 


0.3134 

0.1938 

-0.1451 

-0.4086 

-0.4461 

-0.2730 

0.1938 

0.1641 

0.0121 

-0.1203 

-0.1575 

-0.1016 

-0.1451 

0.0121 

0.3015 

0.4939 

0.4790 

0.2811 

-0.4086. 

-0.1203 

0.4939 

0.9291 

0.9360 

0.5571 

-0.4461 

-0.1575 

0.4790 

0.9360 

0.9519 

0.5685 

-0.2730 

-0.1016 

0.2811 

0.5571 

0.5685 

0.3399 
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Table 3, Eigenvalues of the Closed-Loop System (Example 1) 


Mode i 

Real Part^ 

Imaginary Part 

1 

-1.0000 

±1.183 

2 

-0.4999 

±1.543 

3 

0 

±3.175 

4 

0 

±3.980 

5 

0 

±4.687 

6 

0 

±5.470 


^ Zeroes are less than 10 in double precision 
arithmetic of the AS/9000N system. 


Example 2; 


The second example has a stiffness matrix of 



and natural frequencies given in Table 4. The six lowest 
modes were damped by shifting the real parts of the ei..?en- 
values to -0.5 for each mode. 
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Table 4. Natural Frequencies in Example. 2 


Mode i 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


0.2233 

0.8885 

1.981 

3.475 

5.339 

7.530 

10.00 

12.69 

15.55 

18.51 
21.49 
24.45 
27.31 
30.00 
32.47 
34.66 

36.52 
38.02 
39.11 
39.78 
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The velocity only feedback gain matrix G has been 


turned out to be; 


where 




= 


<=2 = 


°3 ■ 


-0.1214 

-3.1297 

-3.1297 

-1.8767 

-0.2903 

-4.3778 

-4.3778 

-1.7087 

~0.4622 

5.3003 

5.3003 

-1.5374 


The computer output of the actuator matrix Bq and the damping 
matrix C is given in Appendix E. And the eigenvalues of the 
closed- loop system shown in Table 5 indicate that the six 
lowest modes were damped but the rest of them remained un- 
damped as predicted. 
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Table 5. Eigenvalues of the Closed-loop System (Example 2) 


Mode i Real Part^ Imaginary Part 


1 

-0.5006 

±0.6675 

2 

-0.4984 

±0.7856 

3 

-0.5007 

±1.936 

A 

-0.4988 

±2.319 

5 

-0.5000 

±3.123 

6 

-0.4998 

±4.690 

7 

0 

±10.00 

8 

0 

±12.69 

9 

0 

±15.55 

10 

0 

±18 . 51 

11 

0 

±21.49 

12 

0 

±24.45 

13 

0 

±27.31 

14 

0 

±30.00 

15 

0 

±32.47 

16 

0 

±34.66 

17 

0 

±36.52 

18 

0 

±38.02 

19 

0 

±39.11 

20 

0 

±39.78 


— 1 ? 

^ Zeroes are less than 10 in double precision 
arithmetic of the AS/9000N system. 


CHAPTER V 


OPTIMAL CONTROL OF SELECTED MODES 


5 . 1 Introduction 

Consider the state equation, defined in (2,14), for a 

plant 


n(t) = An(t) + Bf(t) 


(5.1) 


where AeR^^^^^^ is the system matrix, is the actuator 

matrix, n(t)eR^’^^^ is the state vector and f(t)eR®^^ is the 
control force vector. Let J be the associated scalar cost 
functional with 


J(n,f,t) = %n (“)Hn(“)+% 


[n^(t)Q^n(t)+f^(t)Q2f(t)]dt. 

(5.2) 


The Hamiltonian for the system is 


H(n,f,r,t) = %n^(t)Q^n(t)+%f^(t)Q 2 f(t)+r^(t) [Ari(t)+Bf(t)] 

(5.3) 


from which it follows that ri(t), r(t) and f(t) must satisfy 
the equations : 


n (t) 

= An(t) +Bf(t) 

(5.1) 

r(t) 

= -Qj^n(t) - A^r(t) 

(5.4) 

0 

= Q 2 f (t) + B^r(t) . 

(5.5) 
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The desired control for minimizing the cost function is 

[10, P.211] 

f(t) = -Q‘^B^r(t) (5.6) 

where it will be assumed that 

r(t) = Pn(t). (5.7) 

Differentiating (5.7) with respect to t and using (5.1) 
and (5.4) we obtain the algebraic Riccati equation, 

+ A^P + PA - PBQ2^B^P = 0 (5.8) 

for P. Substituting P into (5.7) and the resulting equation 
into (5.6) gives the control 

f(t) == -Q"^B^Pn(t). (5.9) 

This control will give the closed-loop matrix 

A = A - BQ2^B^P. (5.10) 

It is usually assumed that H and are sjmimetric posi- 
tive semi-definite matrices and Q 2 is S 3 mnnetric positive 
definite. The Riccati matrix obtained from (5.8) will also 
be s 3 nnmetric and positive definite. The matrices H, and 
Q 2 are weighting matrices chosen to fix the cost penalty for 
the initial conditions, the displacements and the control 
efforts, respectively. 
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This optimal control procedure works quite well when the 
number of modes in the system is not large but the computa- 
tional load for several hundred modes makes this type of ' 
control impractical if time varying gain is used. Even when 
a constant gain is used, the computation of gain is not a 
trivial task. 

There have been niimerous papers about applying optimal 
control theory to the large space structure with the devel- 
opment based on reduced-order models, see [9“, 12 , 13 , 31 ] . The 
computational load can be reduced significantly by this 
approach but the reduced-order model must be carefully cho- 
sen if mode spillover is to be avoided. 

The work in this chapter will take an entirely different 
direction. The computational load for the procedure is 
reasonable and the mode spillover problem can be eliminated. 
The spectral factorization algorithm will be used to decouple 
the selected modes from other modes of the structure. The 
optimal control theory will then be used to construct the 
feedback for the selected modes. The uncontrolled modes are 
uncoupled from the control modes and the possibility of mode 
spillover is eliminated. 

5 . 2 Mode Decoupling of the State Matrix 

In this section a method is presented that decouples 
some of the modes from the remaining ones so that the optimal 
control strategy can be carried out on a lower-order system. 
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This is an extension of the method given in Section 4.1, which 
decouples the stiffness matrix K. 


Consider the undamped system matrix A as defined in 
(2.14) with C = 0 such that 



(5.11) 


where is positive definite. The eigenvalues of K are 

given by where ±jw . are eigenvalues of A. This suggests 
that the spectral decomposition of A can be obtained from 
considering K rather than A. 

Suppose, as shown in Theorem 4.1, that there exists a 
similarity transformation T^ (denoted by T in Theorem 4.1) 
such that 


T^KT 


-1 

k 


^1 ^qx(n-q) 

°(n-q)xq ^2 


(5.12) 


where has q eigenvalues and ti3.s 

(n-q) eigenvalues of K. If such a matrix exists, then Kg^^ 

and Kg 2 give the spectral decomposition of K. To find T^, 

eigenvectors of K must be found (see Lemma 4.2) or the sign 

algorithm [57] can be used to generate T^. The eigenvector 

procedure will probably be the most efficient for large 

systems so the procedure will be based on the method given 



in Lemma 4.2, 


original 
OF POOR QUALfTf 


Let $ denote the orthogonalized eigenvector matrix of 
K and be partitioned with appropriate dimensions as in 
Lemma 4.2, 

*11 *12 
* " *21 *22 


If we further denote orthogonalized matrices of by 
i=l,2, then according to the lemma, K can be block diago- 
nalized to (5.12) by a similarity transformation 

ix 

where 


\ll 

T 

•^K12 

4 

*11*11 

hl*21 


’^K22 


_^22*12 

^22^2^ 


(5.13) 


and 



I . .0.$^. 
11 1 11 


i= 1,2. 


Thus far, the spectral decomposition of K has been car- 
ried out but the system matrix (5.11) must be considered 
because this is the matrix of concern. Let be a new 
transformation matrix with 



(5.14) 
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then, 




'n 


n 




'n 


(5.15) 


This gives a new system matrix with Kg in the lower left 

corner of the matrix, but the matrix (5.15) is not block 

-1 

diagonalized. To block diagonali^.e T^AT^ , construct a 
row- column interchange matrix F, where 


F 


■q I 


L _ I ' 

' In a ~ ~ 

L i _ J 


I 



The blocks of zeroes are deleted in the matrix F for sim- 
plicity. The block diagonal form Ag can then be found by 



TAT'^ , 

where T = FT^ 

and 



0 I * 

q q 1 



-■'bI °q 1 

°2q 

Ag 

r 

0 I 


°2n-2q j 

n-q n-q 


-K.,„ 0 


( 

B2 n-q 


(5.16) 
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OF POOR QUALITY 


It can be shown from (5.13) that T is orthogonal and, after 
simple arithmetic, T is written as 


T 


'^Kll \l2 
0 0 

^21 "^K22 

0 

q n-q 


0 

“^Kll 

0 

%21 

q 


0 

\l2 

0 

^K22_ 

n^q 


}q 

>q 


>n-q 

}n-q 


(5.17) 


where Tj^^j > i,j=l,2 is block partitioned matrices of as 
defined in (5.13). 

The spectral decomposition process will modify the 
state vectot n(t) as shown below. Remember that the state 
vector n(t) was defined in chapter IV as 


n(t) 


x(t) 

x(t) 


(4.14) 


where x(t) is the generalized displacement vector, 
be defined as the transformed vector, 

q(t) = Tn(t) , 

then 

q(t)' = TAn(t) 

but 

n(t) = T"^q(t) 


Let q(t) 


(5.18) 

(5.19) 


thus 

q(t) = TAT"^q(t) . 


(5.20) 
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Therefore, the similarity transformation on A to block diag- 
onalize A will map n(t) into a new vector q(t) as defined 
in (5 . 18) . 

All of the computations for the decomposition given in 
this section can be carried out by considering the K matrix 
which is (nxn) . It is not necessary to find eigenvectors 
of A since the necessary information is contained in. K. 


5.3. Optimal Control of Undamped Decoupled Systems 

It was shown in the previous section that the state 
matrix could be block diagonalized with selected eigenvalues 
of AeR placed in one of the selected block matrices. 

Let the block matrix for the undamped system have the general 
form: 



^B1 

0 


0 }2q 

Ajj 2 }2n-2q 


2q 2n-2q 


(5.21) 


where has eigenvalues |X£l<p and Ag 2 e^^ ^ 

has eigenvalues |X^|>p with p a scalar variable and A is the 
undamped matrix. The value of p will be chosen to include 
the desired modes in A^^. 

Consider now the algebraic Riccati equation for PeR 
and let Ag be the decoupled matrix, thus P must satisfy 
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( 5 . 22 ) 


+ A| f + PAg - PBQ2^B^P. = 0 

where and Q 2 eR^”^ are weighting matrices for q(t) 
and f(t). The matrix BeR^^™ represents the control input 
matrix where 


q(t) = Agq(t) + Bf(t) 


(5.23) 


with B = TB. It is assumed that the algebraic Riccati 
equation (5.22) is completely decoupled such that 


^11 ^Bl^l 
^12 ^ B 2^2 


■f P]_Agj_ 

- PiSiQj’-Sjfi = 0 

(5.24) 

^ 2 ^B 2 


(5.25) 


where Q ^2 the weighting matrices for states and 

= [bJ , 62 ] with and B 2 eR^^^”^^^^. The Riccati 

matrices P^^ and ?2 can be found independently since the 
equations are decoupled. 


Substituting (5.21) for Ag in (5.22), 

+ (T"^)^A^T^P + fTAT"^ - PBQ^^^P = 0 (5.26) 

and rearranging (5.26) gives 

T^Qj^T + A^T^PT + T^PTA - T^PTBQ^^B^T^PT = 0 (5.27) 

Defining P = T^PT and = T^Q^T also gives 
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Q;L + A^P + PA - PBQ'^B^P = 0 , 


(5.28) 


which is a usual algebraic Riccati equation for the general 
control problem. 

Denoting equation (5.24) as system 1 and (5.25) as 
system 2, it follows that system 1 has the system equation 


qi(t) = Agj^qj^(t) + Bj^f(t) 


(5.29) 


where q(t) = 


qi<t) 

q2(i) 


and with cost functional, 


J^(q^,f,t) = 


[qi(t)Qiiqi(t) + f^(t)Q2f(t)]dt, 


(5.30) 


The other system has the system equation 


q2<t) = + B2f(t) 


(5.31) 


and cost functional 





(t)Qi2q2(t) + f^(t)Q2f(t)]dt, 


(5.32) 


where the final state cost was neglected. Assiiming 
that the first system is the desired system for damping, 
then §2 = 0 will leave system 2 undamped and P 2 = 0^ . 

It then follows that the Riccati equation for the uncoupled 
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system is 


ORIGINAL PAGE fS 
OF POOR QUALITY 


P « t^PT <= 


rpL •g m rpt g rp 

nrini •^11^1^12 


ppL g fp mt g m 

n2^in2 


where 


T = 


^11 ^12 


T21 T22 


(5.33) 


and • Since is symmetiic 
P is also symmetric as desired. The control input matrix B 
will have the form 


B 



1 

(2n-2q)xm 


Therefore, since T is orthogonal 


B = 



= T'^B 




L^nxmJ 


(5.34) 


(5.35) 


where B 2 eR’^“. 

^ p 2 is not necessarily identical to zero, but certainly 
§ 2 = 0 is a solution to (5.25) when §2 = 0 and Q ^2 = 0 as 
assumed. This assumption is necessary to obtain a simple 
feedback control law. 



The ntanerical value of can be c’ en so that system 1 is 
controllable and that = 0 but B 2 ^ if the form (5.11) 
is to be maintained. 

The closed-loop system matrix of system 1 after optimal 
control is then written as 

^1 " -^Bl " ®1^2^^1^1 ' 

and, from the equations (5.10) , (5. 11) , (5.33) and (5.35) with 
Bj^=0, the closed-loop system matrix of system (5.1) with 
cost (5.2) and'H=0 is given by 

0 I 

n,, n 

t 

2K+B2Q2^B2tJ2PiT^]_ ®2^2^®2^12^i’^12 

(5.37) 

The required feedback control vector f(t) was obtained 
from (5. 9), (5. 23) and (5.33): 

f(t) = -Q“^B^PTn(t) (5.38) 

which can now be determine J as Q 2 , S and P are known. 

1_ t t“ 

It may be po.v.v’,ib''..i :o make B 2 Q 2 ^2"^12^ 1*^11 “ ^ (5.37) 

by properly selecting B,, and the weighting matrices. In 
general., it will net h<-. ?:•' and the stiffness of the struc- 
ture will be ('hangC'.. . Tc should also be pointed out that 
the matrix C(=B 2 Q 2 ^^ 2 ^ ' ;; doe..- rot r^oresent a model 


A = 
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with passive damping as C does not have the proper structure. 
If the closed- loop system has the feedback defined as (5.38), 
however, there are no restrictions since this control law 
is not for a passive system. 


5.4. An Illustrative Example 
Example 3: 

The example 1 in Section 4.4 will be used to illustrate 
the computational procedure. The stiffness matrix was 


9 -5 

-5 11 -6 

-6 13 -7 

-7 15 -8 

-8 17 -9 

-9 19 


and the mass matrix was Ig. The system was defined in state 

t t * t 

space form with state vector ri (t) = [x (t), x (t) ] : 


n(t) 



»6 ^6 


®1 



n(t) + 



_-K 0g_ 




f (t) . 


(5.39) 


After an equivalence transformation T, which also decoupled 
the two lowest modes, was applied to (5.39), the system was 
changed to 
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q(t) = 


02 I2 


'8x4 


'4x8 


O4 I, 

‘^2 O4 


q(t) + 


B, 


B, 


L 


f (t) , 


(5.40) 


where q(t) = Tn(t) with T given in Appendix F and 


K 


B1 


Kg2 


r 3.83448 

-1.65133 



L^l. 65133 

2.35772_ 



16.81055 

0.23648 

6.93450 

2.07957 

0.23648 

17.86965 

-1.02621 

3.41672 

6,93450 

-1,02621 

24.26231 

-5.16800 

2.07957 

3.41672 

-5.16800 

18.86528 


Thus, system 1 was decoupled and given by 


qi(t) 


“2 12 

■Si s 


q^(t) + B^f(t) 


(5.41) 


which was to be controlled to minimize the cost 


d 1 ( q^^ , f , t ) 


Jqi(t)Q]_j_q3^(t) + f^t)Q2f (t) ] dt (5.42) 


with given “ ^4’ ^2 ~ ^1 B^ = [0 0 -1 1]. The input 
matrix B^ was chosen so that system 1 was controllable. 

The Riccati equation (5.24) corresponding to (5.41) 
and (5.42) was solved by the eigenvector method [61] and a 
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solution was given by 


P 


1 


6.31043 

1.42106 

-0.47965 

-1.24103 


1.42106 

5.51307 

0.90973 

0.84844 


-0.47965 

0.90973 

2.66756 

2.46582 


-1.24103 

0.84844 

2.46582 

4.10804 


With the assumption, B 2 = 0 the solution P of (5.28) 


was computed by (5.33) and listed in Appendix F and by 


(5.36) the closed-loop system matrix of system 1 was 



0 

0 

-4.5959 

2.4127 


0 

0 

1.5900 

-2.2964 


1 

0 

-0.2017 

0.2017 


0 

1 

1.6422 

-1.6422 


Finally, the feedback control vector f(t) was computed 
from (5.38) : 

f(t) = rn(t) 

where 

r = r-0.4845 -0.5051 -0.2486 0.0167 0.1396 0.1097 

-0.1796 0.1450 0.6793 1.0094 0.9489 0.5501]. 


The closed-loop system matrix of the original system was 
given in Appendix F and the eigenvalues of the system before 
and after optimal control were listed in Table 6. 

The input matrix B and the state cost matrix of the 
original system was computed by (5.35) and = T^Q^T. B 2 
is given below and in Appendix F. 
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B 


2 


-0.6696 

-0.4798 

0.1589 

0.6767 

0.7775 

0.4835 


Table 6. Eigenvalues of the System Matrix (Example 3) 


Mode i 

before control 

after control 

1 

±j 1.135 

-0.17l3±jl.l41 

2 

±j2.215 

-0.7506±j2.117 

3 

±j3.l75 

±j3.175 

4 

±j3.980 

±j3.980 

5 

±j4.687 

±j4.687 

6 

±j5.470 

±j5.470 
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CHAPTER VI 


CONCLUSION 


6 . 1 Conclusion 

For the vibration suppression of the large flexible 
space structure, two independent velocity only feedback 
control schemes — the eigenvalue relocation and the optimal 
control — of the second-order system were developed in this 
work. These methods were based on the properties of the 
lambda-matrices and on an efficient mode decoupling technique 
by which selected modes were damped with the rest of the 
modes retaining their pole locations. 

The eigenvalue relocation method allowed the designer to 
place the closed-loop system eigenvalues within the feasible 
region as illustrated in Fig. 1. This development was made 
possible by the aid of the properties of the lambda-matrices 
discovered in Section 3.2. Theorem 3.4 is considered funda- 
mental in this class of system and Lemma 3.5 provides the 
rule on which the eigenvalues move, 

As a result, the vibration control by the eigenvalue 
relocation was accomplished by the following sequence. First, 
through the decoupling procedure described in Section 4.1 
modes were selected for which damping was required. To each * 
decoupled subsystem a damping matrix was computed in such a 
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way that the subsystem had the pre-asaigned eigenvalues. Next, 
under the assumption that the actuators and the sensors are 
collocated, a time -invariant velocity only feedback gain was 
determined in the process which brought the coordinates of 
the subsystems back to the original generalized coordinate 
system. The whole procedure was demonstrated by numerical 
examples . 

Similar results as those achieved in this work may be 
obtained through the pole assignment by gain output feedback 
methods reported in [26-29], but our approach is completely 
different from theirs: our method is devised by using lambda- 
matrix with the constraint that K is invariant rather than 
on the state space which resulted from the conversion of the 
second-order system equations. In addition, it is based on 
the assiimption that the actuator and the sensor matrix are 
designed at our disposal instead of being given as a part of 
a plant. Therefore, it may be meaningless to compare these 
methods and no attempt was made to this end. Roughly 
speaking, however, their methods can assign min(n, m+r-1) 
poles arbitrarily close to min(n, m+r-1) specified symmetric 
values but nothing is said about the remaining poles, where 
n, m and r are the number of states, the rank of actuator 
matrix and the rank of sensor matrix, respectively. Wtiereas , 
in this work 2m poles (m==r and n>2m) can be assigned almost 
arbitrarily within the feasible region as conjugate pairs 
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with the rest of the poles unchanged. 

One advantage of this method over the classical normal 
mode (complete decoupling) technique is that the feasible 
region is much wider than that of the classical normal mode 
case as shown in Figs. 3 and 4. This freedom of choosing 
the locations of system eigenvalues can enhance the damping 
significantly, especially at the low frequency modes where 
more damping is required. 

The eigenvalue assignment technique developed in this 
work also has some shortcomings. First of all, the system 
of simultaneous equations, (4.12), is non-linear. Conse- 
quently, the existence of a solution is not guaranteed, nor 
is Lhe uniqueness of a solution. Nevertheless, a solution 
never failed to exist during the course of .computational 
experiments. Since the damping assignment is achieved at a 
(2x2) subsystem level, computational procedure is not very 
involved, but the simultaneous equations require some 
iterations . 

Another disadvantage of the velocity only feedback 
control scheme is that the sensor and the actuator matrix 
must be computed. When they are given as a part of the plant, 
instead, the gain matrix may be approximated by the least 
square method, for instance. If this is the case, the eigen- 
values of the closed- loop system will certainly deviate from 
the pre-determined position even though the error may not be 
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significant because of the continuity of the eigenvalues. 

This problem remains for further research. 

On the other hand, a procedure for optimal control of 
selected modes was also developed in the last chapter 
and it was shown that the control vector for a rather small 
second-order lambda-matrix could be determined in such a way 
that damping was added to the lowest modes with the other 
modes remaining unchanged. The elimination of mode spillover 
problem was made possible by the mode decoupling procedure 
and by the manipulation of the actuator and the sensor 
matrices . 

Finally, the decoupling procedure for both methods 
requires only the eigenvectors of the assigned modes. This 
fact resulted in a significant saving of computing time when 
the number of modes involved in damping assignment was less 
than one fourth of the total modes . 

The main contribution of this work can be summarized as 
follows: i) unique properties of the second-order lambda- 

matrix were discovered and applied to the vibration control 
problem, ii) a computational procedure for the damping 
matrix determination with the stiffness matrix invariant was 
established, iii) a technique to decouple the large system 
into smaller .^iubsystems through partial eigenvectors was 
developed, and iv) an optimal control method of the selected 
modes- without any spillover to other modes was investigated. 
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6 . 2 Recointnendations for Further Work 

One of the topics which deserves further research is the 
passive damping algorithm: the complex sensory and actuating 

devices of the active control may provide enough motivation 
for research on the passive damping of the large space struc- 
ture. In order to realize the passive damping, however, the 
structure of the damping matrix should be in a simple form 
such as diagonal or tridiagonal. Another promising area is 
to extend the theory developed in this work to higher -order 
systems; third-order, fourth-order, and so forth. This ex- 
tension may be established without much difficulty since 
there is well-developed theory on the corresponding lambda- 
matrices. An algorithm of damping matrix determination for 
bigger than a (2x2) subsystem is worth investigating. As to 
the optimal control part of the work, it would be worthwhile 
to find a necessary condition for the Riccati equation to be 
decoupled completely. Finally, it may be interesting to 
combine the eigenvalue relocation technique and the optimal 
control method in such a way that the selected modes are con- 
trolled optimally in some sense with desirable pole locations. 
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APPENDIX A 

LAMBDA- MATRICES AND GENERALIZED LATENT VECTORS [51] 

Let A^eC^^^, i == 0,1,..., m and XeC. A X-matrix is 
defined as 

A<X) = AqX™ + A^x““^ 


A latent root is also defined by a scalar X^eC such that 
ACX^) is singular. When Aq is non-singular the X-matrix 
is called regular and in such a case a monic X-raatrix, • A(X) 
can be obtained by changing coordinates, i.e., 

I 

A(X) = A“^A(X) 

= IX® + A^X®"^ +...+A^ . (A. 2) 

Associated with A(X) a block companion matrix is 
defined as - 



• • 



• « 



(A. 3) 
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There are very useful relationships between eigenvalues and 
eigenvectors of the block companion matrix and latent roots 
and latent vectors, which will be presented in the following 
definition and theorems , 


Theorem A.l . Let A(X) be defined as in (A. 2) and A^ as 
given in (A. 3) then the latent roots of A(X) are the eigen- 
values of A . 

c 

Proof: When (IX - A ) is post-multiplied by the following 

unimodular block Toeplitz matrix. 


V(X) 


n 

XI 


n 




n 


0 


n 


n 


n 

XI 


n 


n 
I 

n n 


the following result is obtained 


(IX-A^) V(X) = 


n 


n 


-I 


n 


n 


0 • • 

n 

A(X) * 


n 


-I 


n 


n 


n 


n 


-I 
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in which A(X) is as given in (A. 2). Now when the determinant 
of the above expression is taken by expanding minors along 
the ones, "it follows that 

det[(lX-A^)V(X)] = det(IX-A^) = (-1)°^ det A(X) 


which implies that the eigenvalues of A^ are latent roots 
of A(X) . VVV 


Definition A.l . Let A(X) be defined as in (A. 2) with 
latent roots X^ of multiplicity n.. Primary right and 
left latent vectors, eC , respectively, are 

defind as 


A(X.)yp> 


= 0 T 

nxl 

= Vi 


j ~ 1 , 2 , > . . , 
j ~ 1 > 2 , . . . , q^ 


(A,4) 
(A, 5) 


where the number of primary right or left latent vectors of 
X^ denoted by q^, is the nullity of A(X^) 

It is known [56] that a primary right latent vector. 
yP^ is a subvector of a right eigenvector of the block 
comp-^.nion matrix (A. 3) with 


y(j) = 

^Ci 


y<j> 

X.y^) 




(A. 6) 
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where ^ satisfies the algebraic equation (A - = 

^mnxl* also not difficult to show that a primary 

left latent vector, satisfies a similar form as the 

primary right latent vector with 


(j) ^ 


cx 


a'l-h + +...+A‘_2)zp> 




z 0 ) 

i 


(A.7) 


where is an eigenvector of the block companion matrix. 

It should be noted that the maximum n^Jmber of primary 

right or left latent vectors is n, i.e., max(q.) == n. If 

i ^ 

n. > q. then n. - q, generalized latent vectors must be 

X. JL Jiv 

constructed to define the complete set of latent vectors of 
the lambda matrix. These generalized latent vectors can be 
obtained from a chain rule given in [57]. Notice that each 
of these n^ - q^ generalized latent vectors may or may not 
have a chain of vectors according to the structure of the 
Jordan block of the block companion matrix. The length of 
the chain for each primary latent vector could be determined 
during the computation of vectors until no vector satisfies 
the chain rule. 
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Theorem A. 2. Let A(X) be defined as in (A. 2) and a set 


of right latent vectors form a 

right Jordan chain associated with the latent root and 
the j-th primary latent vector. Then, the chain rule is 
given as 


A(X )yt« + + 1 ^ y<<t-2) 


+ 


(A-1) I , (£~r) ^i " ^ 


nxl' ^ 1,2, . , . ,hj 

(A. 8) 


where y^^^ is the j-th primary right latent vector and h. 

■L J 

is the length of the Jordan chain. The vectors y^ for 
1 < k < h are generalized right latent vectors of the j-th 

- j 

primary latent vector. 


Proof. The proof of this theorem is obtained from consid- 
eration of the chain rule for generalized eigenvectors. 

The chain rule is 


(A^-x^Dy 

(A^-X^I)y 


(A -X.Dy 
c 1 


(2) 

- 

ci 

y ci 

(3) 

- v<2) 

ci 

” ^ci 

• 

(h.) 

« 

eiJ 

- 

- y<.iJ 


(A. 9) 
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being the j-th primary eigenvector and y^^^ 'for 

1 < k < h. its generalized eigenvectors associated with the 

bJ 

eigenvalue of A^. When (A. 9) is expanded the chain rule 

/Ir) 

is obtained where y^ is formed from the first m rows of 



Generalized right latent vectors yf^^ and will 

i 1 

be obtained from 


(2) = 0 


A(J )y(3> + y<-2> + i ^ „(1) 

AU^Jy^ ^ . dx ^i ^21 ^,2 ^i 


= 0 


dx‘ 


nxl 


or (A. 8) in general. VVV 

The chain rule can also be utilized to modify (A. 6) for 
the relation between generalized eigenvectors and general- 
ized latent vectors with 



.(k) 




x-' X X*' X X 


m-1 

I 

j=0 


( 


n-1 

j 




(k-j) 

^i 


(A. 10) 
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Theorem A, 3 . Let A(X) be defined as in (A. 2), a set of 
left latent vectors » • • • from a left 

Jordan chain associated with the latent root and the j -th 
primary latent vector. The chain rule is given by 


2*t, 


tp I ^ ^^i^_(A-i) . 1 ^'^i^ _a-2) . 

(X,.)z,. + ^ Z^. + -^ Z,. +... 


1' i 


2' ^ 


■*■ ( 1 - 1 ) 1 


z'l) 

L 


= 0 


nxl 


, 5- = 1 , 2. . . ,h. 


(A. 11) 


where z^^^ is the j-th primary left latent vector and h. is 
•I* J 

/V) 

the length of the Jordan chain. The vectors z) for 
1 < k < hj are generalized left latent vectors of the j-th 
primary Is tent vector. 


Proof; The proof of this theorem follows directly from 

fk) 

generalized left eigenvectors z^^ of the companion form. 
Generalized left latent vector z^^^ will be formed from the 


last m rows of . 

cx 


VVV 


,(k) 


Generalized left eigenvector z^^'' can also be defined 

from latent vectors z^^^ of A . , the latent roots A. and 

X 1 i 

the lambda matrix A(A) . Utilizing (A. 7) and the chain rule 

(k) 

for z^^ , it follows that left generalized eigenvectors 
satisfy the relation; 
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^ci 


\iii— 1 T , •t»iH“2 


HAjX“-‘+. . .4A;.j)zf >+(m-l)X‘“-"l+(m-2)A^Af^+. . -•'4. 


m-2. 




,(k“l). 


(X^I-HaJx^+A^) z^^4(2X^I+aJ) 


Ak) 


(A. 12) 


where the latent vector 



is defined only for 1 < k < h. . 
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APPENDIX B 


SOLVENTS OF MATRIX POLYNOMIALS 

Associated with an m-th order lamb da -matrix (A. 2) two 
types of pol 3 niomial, the right matrix pol 3 niomial 


Aj^(X) = X”^ + A^X“ . .+A^, 


(B.l) 


and the left matrix polynomial 


Aj^(X) “ X°^ + x”"^Aj+. . .+A^ 


(B.2) 


are defined for XeC’^^’^. Matrices Xj^, X^eC^^^ called a 
right and a left solvent if = 0^^^ and Al<Xl) = O^xn* 

respectively. In this appendix structures of solvents will 
be examined in terms of latent roots and latent vectors and 
theorems on the existence of solvents will be given at the 
end. 

n.xn. 

Theorem B.l . Let J^eC ^ be a single Jordan block with 

latent root X. of multiplicity n. < n and y^^\..., 

(ti.) nxl “ ’ i 1 

y^ eC be generalized right latent vectors of A (A) , then 


m 


k=0 ^ ^ ^nxn. 

3 . 

,( 1 ) „( 2 ) 


(B.3) 


(n. ) 


where ‘ ' ^i ^ ^ ^0 
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/IN S J \ 

Proof; Let = ^Yci^ ’•' ^ci where is a primary 

right eigenvector of the block companion matrix A (A. 3). 

c 

Then it is not difficult to obtain Y . from (A. 6) and Y . 

^ ^ iL 

is written as 


Y , 

Cl 




Y. 

i 1 


(B.4) 


From the chain of generalized eigenvectors (A. 9) of the 
block companion matrix A^ it follows tha;; 


i.e 


f 


C Cl 


Y .J. 

Cl i 


0 _ I_ f', , • • • • 0 ,„ 


Y. 

n n :i n 


1 

[ — 1 
o 
o 


Y.J. 

n n - n 

« » • • 

• • • * 

B • • • 


i 1 

• 

* 

’^1 


11 

_ 





Vi 

Y 


Y.j“ 
1 1 


(B.f) 


Thus, the last row of (B.5) provides the conclusion of the 
theorem, VVV 


n .xn , 


Theorem B . 2 . Let J^eC ^ ^ be a single Jordan block with 

latent root X. of multiplicity n.< n and 
(n.) ^ ^ ^ 1 - 1 


( 1 ) 


> • • • I 


- - w *1 

3^ ^ be generalized left latent vectors of A(X) , then 
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n 


V T^i-ky 

-n.xn 


,A,. = 0„ 


k=0 


(B.6) 


(n. ) (n. -1) ^ 

where A« = I and Z. = [ 2 . ^ z ^ ... 2 : M . 

V/ ^ JL ^ M 

(n.) (n.-l) Qv ^ 

Proof*. Let Z . = [z . 2 . .... , where z^i'' is a 

0«!« ClCiL C*l> C>L 

primary left eigenvector of the block companion matrix 

A^(A.3). Then it follows from (A. 12) that 


Z^l = [2. J.2^ ... Jfl 2.1 


(B.7) 


and 


Z^.A = J.Z^. 
Cl C 1 Cl 


1 . e. 


[Z. J.Z j“"^Z.] 

Ill 1 


o„ 

0„. 

• 0„ 

-A 

n 

n 

n 

m 

I 

0 . 

.0 

-A 1 

m-1 

n 

n 

• 

« 

* 

n 


On- 

•In 

1 

1 


= [J.Z. J?Z J?"^Z. ] 

1111 1 1 


(B.8) 


The proof is obtained from the 





n.j^xni 

Theorem 

B.3. 

Let 

J.eC b( 

1 

in Theorem B . 1 , 

Y = 
SL 

[Y^ Y 2 ... 


with 

i 

E n. 

=1 1 

= n . Then 


On 

In 

On • • • 


On 

On 

In • * • 


~A “A 1 "A 9 • 
m m-i m ~2 


last column of (B.8). VVV 

! defined as the same way as 
Yj] eC"™, and J = 


~ 





On 


Y 


YJ 

On 


YJ 

« 


2 

YJ 

A, 


• 

4 

Yjm-1 


Yj“ 

Li 

- 

L J 


iT' a. 
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Proof: Since 


~ n 




Y 




YJ 

« ■ 


•flJl 

• 

» 

r 

s 

1 

1 


« 

* 

Y J^"^ 
jri 

• 

Y’jin-i 

’ • ■ t, 


the proof is completed when (B.5) is applied block column- 
wise A times. VVV 

n.xn. nxn. 

Theorem B . 4 . Let J^eC , Z^eC ^ be defined as the 

same way in Theorem (B.2) and Z = [2^ Z^ ••• 2^]^, 

J => Bdiag(Jj^, J 2 , ^ . ,dj^) with^^^n^ = n, then, 


[Z JZ 


,J°^‘^Z] 


0„ • • • 

0„ 

-A 

n 

n 

m 

I ' ‘ • 

o„ 

-A T 

n 

n 

m-1 


0 ■ ' • I 
\i n 


“A-, 


[JZ J“Z ... J®Z] 


(B.IO) 


Proof; Since 


[Z JZ . . .J°^'^Z] 


J.Z, . . . J?”^Z, 
1 1 11 

T 7 ^ 7 

J 2^2 • " 2 2 


7* J *7 J°^~^Z 
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when (B.8) is applied block rowwise times the proof is done. 

VVV 

Theorem B . 5 . Let be the same as in Theorem B.3. 

If n columns of (n right latent vectors including 

generalized right latent vectors) are linearly independent, 
then YJY”^ is a right solvent of A(A). 

Proof: From the last block row of (B.5) 

m . 

with Ao = 1. (B.ll) 

Since Y "*■ exists by the hypothesis post-multiplying y"^ at 
both sides of (B.ll) gives A^(YJY ^) = which implies 

that YJY~^ is a right solvent of A(X) by the definition. 

VVV 

Theorem B , 6 . Let be the same as in Theorem B.4. 

If n rows of (n left latent vectors including gener- 

alized left latent vectors) are linearly independent, then 
Z~^JZ is a left solvent of A(X). 

Proof: From the last block column of (B.IO) 

Jo withAo = I. (B,12) 

-1 -1 
Since Z " exists by the hypothesis pre-multiplying Z at 

both sides of (B,12) gives A_ (Z“^JZ) = > which implies 

that Z”^JZ is a left solvent of A(A) by the definition, 
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APPENDIX C 
PROOF OF THEOREM 3.3 

Proof [A2, PP. 48-49] We have 

A(X) = A.A^ + A,X’^''^+. . .+A 
u i m 

= A^X^'^(IA-R) + (AqR+A^)x“"^ + A2X““^ 

= [AqX™"^ -F (AqR+A3^)X™"^] (IX-R) 

+ (AqR^+Aj^R-FA2)X®“^ + 

= [AqX”"^+(AqR+A^)X°^"^ +. . . + (AqR™"^+A^R^'^+. . ,+Ajj^_^)](IX 
+ (AqR®+A^r“"^ 

If R is a right solvent of A(X) 

A^(R) = A.R”^ + A,R°^‘^+. . .+A = 0 

K 0 1 m nxn 

Thus, the last equality of (C.l) implies A(X) is divisible 
on the right by (IX-R) with quotient of (m-l)th order 
lambda-matrix. This provides the necessity of the theorem. 

On the other hand, if A(X) is divisible on the right 
by (IX-R) Aj^(R) must vanish, which is the definition of a 
right solvent. Now, the sufficiency of the theorem is proven. 

In a similar manner, divisibility on the left by (IX-L) 
can be shown without any difficulty. VVV 
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APPENDIX D 

PROOF OF THEOREM 3.11 [59] 


For an arbitrary matrix A, define 

i |(A + A*) 
i ^(A - A*) 

then, H 2 are Hermitian and 
A = + jH2. 


According to the Rayleigh Principle, for every eigenvalue 
X of A we have ? 


ReX 


max 

x=s^0 


Re 


x*A* x 

x*x 


1 1 

max -^(x*Ax + x*A*x) 
x?^0 ^ ^ ^ 


TT13X 

x’^O 


x*Hjj^x 

x*x 


^max 


(Hj^) . 


Similarly, 


and 


ImA 

ImA 

ImA 


< 

> 


max 

x^O 

min 

x=/0 



x*Ax 

x*x. 


Re 


x*Ax 

x*x 




Vin^^l^ ' 


> 

- mxn 
x^^O 


Im 


x*Ax 

x*x 


mxn Z 


VVV 
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APPENDIX E 

THE ACTUATOR MATRIX Al^D THE DAMPING MATRIX C OF EXAMPLE 2 


ACTIIATilR matrix 


-0.aS99SI)-01 

-0.173flflD+0n 

0,909630-01 

’',209910+00 

•0,133900+00 

0.241260+00 

-0. 909630-01 

-0,28727r+00 

0.17364O+U0 

0.307740+00 

-0,241280+00 

0,300670+00 

-0.133900+00 

-0.300670+00 

0.241260+00 

0,241280+00 

•0.300870+00 

0.133900+00 

-o.!73aap+oo 

-0, 209910+00 

0,287270+00 

0,459950-01 

-0.300870+00 

-0, 133900+00 

-0.20991D+00 

-0.45995P-01 

0,307740+00 

•0.173840+00 

-0,241280+00 

-0.300870+00 

-0,201280+00 

0.13390D+00 

0,300870+00 

-0.300870+00 

-0,133900+00 

-0.241280+00 

-0,267260+00 

0.267260+00 

0.267260+00 

-0.267260+00 

0,476330-14 

0.409230-15 

-0.287270+00 

0.307740+00 

0.209910+00 

-0,909630-01 

0.133900+00 

0.241280+00 

-0,30087D< jo 

0.2412RO+00 

0,133900+00 

0.133900+00 

0,241280+00 

0. 300870+00 

-0.30774t-+00 

0.909630-01 

O.4S99SD-01 

0. 287270+00 

0.300870+00 

0.133900+00 

-0.307740+00 

-0,909630-01 

-0,459950-01 

0,287270+00 

0,300870+00 

-0,133900+00 

-0,300670+00 

-0.241280+00 

-0.133900+00 

0.133900+00 

0.241280+00 

-0.300870+00 

-0.297270+00 

-0.307740+00 

-0. 209910+00 

-0.909630-01 

0,133900+00 

-0,241260+00 

-0. 267260+00 

-0.267260+00 

-0.267260+00 

-0. 267260+00 

-0.547480-14 

0.077070-15 

-0.201280+00 

•0.133900+00 

-0. 300870+00 

-0.300670+00 

-0,133900+00 

0,241260+00 

-0.20991D+00 

0,459950-01 

-0.307740+00 

-0,173840+00 

•0,241280+00 

0.300870+00 

-0.173800+00 

0.209910+00 

-0.287270+00 

0.459950-01 

-0.300870+00 

0.133900+00 

-0,133900+00 

0.300670+00 

-0.241280+00 

0. 241260+00 

-0,300670+00 

-0.133900+00 

-0,909630-01 

0.2BT27D+00 

-0.173840+00 

0.307740+00 

-0,241260+00 

«-0. 300070+00 

-0. 459950-01 

0.173U4O+0O 

-0.909630-01 

0.209910+00 

-0.133900+00 

-0,241280+00 

COMpUTfO OAMPIRS 

MATRIX C ^ 






0,7922 

1.2310 

1.1739 

0.7593 

0.2947 

0,0451 

0.0716 

0,2240 

0,2615 

0,1286 

•0.1622 

-0.3950 

-0,4049 

-0.1824 

0,1237 

0.3199 

0,3082 

0,1454 

-0,0181 

-0.0648 

1.2310 

1.3941 

1,7640 

1,0720 

0,3231 

-0,0491 

0.0527 

0.3760 

0.5468 

0.i589 

•0.0861 

-0,4779 

-0,5441 

-0,2541 

0.1706 

0.4349 

0,3895 

0,1275 

-0,1131 

-0,1501 

1.1739 

1,7640 

1.5468 

0.7784 

0.0007 

-0.3243 

-0.0984 

0.3961 

0.7352 

0.6647 

0,2579 

-0.1752 

-0.340? 

-0.1746 

0.1302 

0.3004 

0.2005 

-0,0709 

-8.2763 

•0,241.9 

0.7593 

1.0720 

0.7784 

0.1004 

-0.5147 

-9,6898 

-0.3505 

0.2559 

0.7503 

0.8780 

0.6505 

0,2909 

0.0434 

-0.0013 

0,0602 

0.0594 

-0.0851 

-0.2954 

-0,4016 

-0.2888 

0,2947 

0.3231 

0.0007 

-0,5147 

-0,9230 

-0.9722 

-0.6135 

-0,0210 

0,5258 

0.8109 

0.7893 

0.5731 

0.3245 

0.1428 

0.0249 

-0.0867 

-0.2222 

-0.3410 

-0.3616 

-0.2354 

0.0451 

-0.0491 

-C.3243 

-0.6898 

-0.9722 

-1.0178 

.0.7821 

-0,3535 

0.0996 

0,4199 

0.5339 

0.4o82 

0,3094 

0,1445 

0.0206 

-0,0578 

.0.1041 

-0.1248 

•0.1148 

-0,0698 

0,0716 

0,0527 

-0,0984 

-0,3505 

-0.6135 

-0,7821 

-0,7870 

-0.6282 

-0.3735 

-0.1232 

0,0382 

0,0807 

0.0365 

-0.0241 

-0.0351 

0.0280 

0.1366 

0.2273 

0.2401 

0.1539 

0,2240 

0,3760 

0.3961 

0,2559 

•0,0210 

•0.3535 

-0.6282 

•0.7505 

•0,6953 

-0,5226 

•0,3424 

-0.2460 

-0,2476 

-0,2761 

-0,2279 

-0,0469 

0.2239 

0.4561 

d.5l23 

0.3370 

0.2815 

0.5468 

0.735? 

0,7503 

0.5258 

0,0996 

-0.3735 

-0.6953 

•0,7448 

-0,^583 

•0.3082 

-6.1651 

•0,2630 

-0,4463 

-0.5395 

-0,3951 

.0.0327 

0.3581 

0.5398 

0,3917 

0.1286 

0.3589 

0.6647 

0,8780 

0.8109 

0,4199 

-0.1232 

•0.5226 

•0,5583 

•0,2485 

0,1507 

0,3195 

0,1018 

-0,3790 

-0,8029 

-0,8800 

-0,5569 

-0.0535 

0.2975 

0,2921 

•0.1622 

-0,0861 

0.2579 

0,6505 

0,7893 

0.5339 

0.C382 

-0.3424 

•0,3082 

0.1507 

0.7247 

0.9892 

0.7020 

-0.0208 

-0,7800 

-1.1643 

— 1 1 1)226 

-0.5427 

-0,0867 

0.0661 

•0,3950 

-0,4779 

-0.175? 

0,2909 

0.5731 

0.4682 

0,0607 

-0.2460 

-0,1851 

0,3195 

0,9892 

1.3891 

1.2228 

0.5335 

-0.3280 

-0,9399 

-1.0753 

-0.8120 

-0,4214 

-f ,1391 

'0,4049 

-0.5441 

-0.340? 

0.0434 

0.3245 

0.3094 

0.0365 

-0,2476 

•0,2630 

0,1018 

0.7020 

1.2228 

1.3755 

1.0712 

0.4596 

-0.1790 

-0.6020 

-0.7178 

•0.5826 

-0.3133 

•0.1824 

-0.2541 

-0,1746 

•0.0013 

0.1428 

0.1445 

-0.0241 

-0.2761 

-0,4463 

-0.3790 

'0.0208 

0.5335 

1.071? 

1.3580 

1.2536 

0,7876 

0.1538 

-0,3777 

•0.5890 

-0,4208 

0.1237 

0.170,6 

0,13p? 

0.0602 

0.0249 

0.0206 

-0.0351 

-0,2279 

•0,5395 

-0,8029 

•0,7800 

-0,3280 

0,4596 

1.2536 

1,6567 

1,4470 

0,7340 

-0.0870 

-0.5694 

-0.4922 

0.3199 

0.4349 

0.3004 

0,0594 

-0.0867 

-0,0578 

0.0280 

-0.0469 

-0.3951 

-0.8800 

'1.1643 

-0.9399 

-0.1790 

0.7876 

1.4470 

1.4321 

0.7748 

•0.0989 

-0,6451 

-0.3641 

0.3082 

0,3695 

0.2O05 

-0.0651 

-0.2222 

-0.1041 

0.1396 

0.2239 

•0,0327 

-0,5569 

■1.0226 

-1,0753 

-0,6020 

0.1538 

0,7340 

0,7748 

0.2663 

-0.4290 

-0.8213 

-0,6343 

0.1454 

0,1275 

-0,0709 

-0,2954 

-0,3410 

-0,1248 

0.2273 

0.4561 

0,3581 

•0.0535 

'0,5427 

-0.8120 

-0,7178 

-0.3777 

-0,0870 

-0.0989 

-0.4290 

-0.8311 

-0,9675 

-0.0536 

'O.oiBl 

-0.1131 

-0,27b3 

-0,4016 

-0,3616 

•0.1148 

0.2401 

0.5123 

0,5396 

0.2975 

0,0867 

-0.4214 

-0,5826 

-0.5890 

-0,5694 

-0.6451 

-0.8213 

-0.9675 

-0,9090 

-0,5595 

0.0648 

-O.ISOI 

-0,2419 

-0.2888 

-0,2354 

-0,0698 

0.1539 

0.3370 

0.3917 

0,2921 

0.0861 

-0.1391 

-'5.3133 

-0.4206 

-0.4922 

«n,'564i 

•0.6343 

-0.6536 

-0.5595 

-0.3273 


Two consecutive lines make a row of the matrix C. 
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APPENDIX F 

MATRICES T, P, S AND OF EXAMPLE 3 

TPAhSFDRMilTIOK matrix T 


0.63876 

0.64973 

0,29040 

-n.070t>7 

•0.22765 

-P.lo93i 

O.OOOOO 

o.oou'jo 


j.*jUOuO 

t'at 


•0.03007 

0.16995 

0 , 4493 a 

0.60598 

0.54966 

0.31416 

t .noooo 

O.uunoO 


1*. IW'J ^ 1 

t • V 

• Uutf Ov 

n.ocooc 

.00000 

O.oonoo 

O.OOOOO 

o.ouooo 

J.oooou 

11.63878 

0 , 04973 


• 'l,t>7i.o7 


16931 

0.00000 

0.00000 

O.nonon 

u.coooo 

O.OOOOO 

O.OuOUl) 

.0,03087 

0.1O4V5 


•1.60898 


v.31916 

•0.08970 

0.07674 

0,7o*an 

-0.49496 

-O.oulis 

-0.07611 

O.nuOUD 

0.00000 

v,«*uuu.) 

•l,OUUo‘* 


v.ObOLO 

0.3S59B 

•0.67761 

0 . 45947 

0.24175 

•0.37430 

•0,06662 

o.noouo 

o.Quoua 


u.oOnoO 


k.bbUUb 

0.41620 

•0.28422 

0.02851 

•D. 43 YI 6 

0.70593 

-0.23757 

o.ouono 

o.uuooo 

il*uUUOO 

O.OUOOO 

i-, 4.*t#gvu 

w.uuiibb 

0.22639 

•0.05608 

•0.00068 

•0.36504 

-0. 08002 

0.89774 

C.noouo 

o.ououc 

U.MUUUO 

O.ouuoi- 

i • ViTuWII 

u.lluouu 

0.00000 

0.00000 

O.oooon 

O.OvOOO 

0.00000 

0.00000 

.0.46970 

0,0767-1 

u*7utS(io 

.0.49490 

— ."HllS 

-u.b/ol 1 

O.OOOOO 

0.00000 

o.nonon 

O.OnOuO 

o.ooono 

O.OOOOO 

0.35593 

-0,677oI 

0,4|b()47 

u.24| 75 

•4 ,3743u 

-•1,40062 

0.00000 

0.00000 

o.oooon 

o.oonoo 

O.OOOOO 

o.ooooa 

0.41620 

-0.26422 

O.D2oSI 

.0.43716 

«',/u543 

-0.23757 

0.00000 

0.00000 1 

o.oooon 

O.OOOOO 

0,00000 

0.00000 

0.22939 

-0,05608 

•U.OUtibA 

•0.365u4 

•I ,u6uu2 

b, 69774 

full RICCAT] 

■ MATRIX P 











2.52410 

2.71587 

1,48925 

0.16516 

-0,50212 

- 0.44334 

.0.18837 

-0,35650 

•0.46512 

•0.4/262 

-u,374ln 

•J.2bUe5 

2.71587 

3.13700 

2.09668 

0.82C46 

0.03452 

-0.15064 

.0.07967 

-0.21456 

-0,34313 

-0.39015 

•0,32634 

•b, 16143 

1.48925 

2.09669 

2.01617 

1.57661 

1.02652 

0.46953 

0.17150 

0.17664 

0.00763 

.0.00641 

• U.-14 969 

•u.b3906 

0.16518 

0.82046 

1 .5761,1 

1.93427 

1.66725 

0,94771 

0.35522 

0.46249 

U. 44037 

0.32335 

0.197/1 

il.uVvOO 

•0.50212 

0.03452 

1.02652 

1.68725 

1.63813 

0.96164 

0.36616 

0.5232b 

0.51355 

0,41064 

0.27314 

u. 13214 

•0.44334 

•0.15064 

0,48953 

0.94771 

0.96160 

0.57383 

0.21972 

0,31946 

0.3207e 

0.2i>291 

V.17055 

u. 06761 

•0.18837 

•0.07987 

0.17isn 

0.35522 

0,36616 

0.21972 

0.99512 

1,30380 

1.12359 

0.70259 

0.42577 

0,17936 

•0.55650 

•0.21456 

0.17864 

0.46249 

0.5232b 

•1.31946 

1.30360 

1.76932 

1.65663 

1.24183 

u, 77467 

U.3So2b 

•0.46512 

•0.34313 

0.08763 

0.44037 

0.51355 

0.32076 

1.I23S0 

1,65663 

1.69793 

1.41947 

0.96015 

0. 46012 

•0.47262 

•0.39015 

•0.00641 

0.32335 

0.41084 

0.26241 

0.76259 

1,24163 

1.41947 

1 ,31065 

0,97575 

0.50625 

•0.37418 

•0.32834 

-0. 04964 

0.19771 

0.27314 

0.17855 

0.42577 

0.77467 

0.96015 

O.V7575 

0.70297 

b. 40654 

•0.20085 

•0.18143 

-0.03908 

0.09000 

0.13214 

0,06761 

0.17938 

0,35626 

0.46613 

0.5U625 

U.4U654 

0.21960 

CLOSEO^LOQP 

SYSTEM matrix t 










0.00000 

0.00000 

0,00000 

O.OOOOO 

o.oonoo 

0,00000 

l.ououn 

O.OUOOO 

O.OOOOO 

C.oupbu 

•u.booua 

•O.OUOOO 

0.00000 

0.00000 

o.oooon 

0.00000 

o.oonoo 

o.oooon 

O.OOOOU 

i.nooou 

U.OoilOO 

-u.cuouo 

u.uubou 

O.UOOOO 

0.00000 

0.00000 

O.oooon 

0.00000 

0.00000 

4,00000 

O.OOOOO 

O.OOOOO 

l.UUOOO 

-n.OuUoo 

b.bOnuo 

0.00000 

0.00000 

0.00000 

o.oooon 

0.00000 

O.OOOOO 

O.OOOOO 

O.OOOOO 

•O.OOOOO 

-O.OUOOO 

l.ouuoo 

b.JOOuo 

-U.OOOOO 

0.00000 

0.00000 

o.oonoo 

0.00000 

0.00000 

0,00000 

.O.ooOoo 

O.OOOOO 

O.OOUOD 

u.nOJuo 

I. 1 . 0 CUD 

•U.OOOOO 

0.00000 

0.00000 

0.00000 

O.OOOOO 

O.OOOOO 

0.00000 

-O.OOOOO 

0,00000 

0,00000 

-u.nuaiiO 

-u.ooobo 

l.bouuu 

>9.32442 

4.66175 

•0.16651 

0.01116 

0.09350 

0.07343 

.0.12025 

0,09912 

0.45492 

0.6/595 

0,63544 

0,3683b 

4.76757 • 

.11.24234 

5.KBQ71 

0.00800 

0.06699 

0.052ol 

.0.08615 

0,07102 

0.32593 

0,46429 

4,45527 

0.26391 

0.07700 

6.08028 -12.96048 

6.99735 

•0.02219 

-0,u1743 

U, 02854 

-0,C2353 

•U. 10798 

.0. |6U44 

-u. 15062 

•U.U6743 

0.32781 

0.34178 

7.16825 

.15.01128 

7.90552 

-0,47420 

0.12151 

•0, 10016 

•0.45967 

.ii.663u2 

-U.042U9 

-0.3/221 

0.37667 

0.39273 

0.1933? 

7.98704 

-17.10856 

8.91474 

0,13962 

-0,11509 

-0,52619 

.0,76463 

-0, 737/9 

-•»,427o9 

0,23422 

0.24420 

0.12021 

•O.OO 8 O 6 

8.93249 

•19.05301 

0.48662 

-0,07156 

-0.32643 

-u.aeoul 

•0,458/7 

•0.26594 

SI CIF 0RI61NAL SYSTEM 











0.40899 

0.40979 

0.17163 

•0.06385 

•0.16239 

-0.11785 

O.OUOOO 

0,00000 

o.ooono 

u.ououo 

n.ouuuu 

u.ouuoo 

0.40979 

0.45103 

0,26505 

0.05707 

- 0.05440 

.0.05661 

O.OOOOO 

0,00000 

o.ouooo 

o.oouro 

o.nvtiuo 

U.UOOOO 

0.17163 

0.26505 

0.28624 

0.25177 

0.18097 

0.09200 

o.ouooo 

o.ouooo 

0. 00(100 

o.ouuon 

0 . .>oOb>) 

u.uOOOO 

•0.06385 

0.05707 

0.25177 

0.37221 

0.34929 

0.20234 

o.noouo 

o.ouooo 

o.uuouo 

U.QUOUO 

u.uoObO 

O.OUOOO 

•0.16239 

•0.05446 

0.18097 

0.34929 

0,35417 

0.21129 

fl.UUOuO 

o.ouooo 

u.oooon 

U.OOOOO 

v.UUUUO 

U.UOOOO 

•0.11785 

•0.05661 

0.092on 

0.20234 

0.21129 

0.12736 

o.noouo 

u.coooo 

o.ouooo 

O.ljUUOO 

L •V.UUUll 

o.iioilOb 

0.00000 

0.00000 

0.00000 

0.00000 

O.OOOOO 

o.oooon 

0.40899 

0. 40979 

u.i7i63 

.•1,00365 

-••, 10239 

-0.11765 

O.OOOOO 

O.OOOOO 

O.OOOOO 

O.OOOOO 

o.ouooo 

Q.OOOOO 

0,40979 

0,45103 

0.2o5ilS 

U..15/v7 

-U,-‘5446 

-u.obaol 

0.00000 

0.00000 

O.oooun 

O.OOOOO 

O.OnOoO 

o.noooo 

W. 17163 

0,265u5 

• 1 . 20024 

0.25177 

•’,leu97 

U.U92UV 

0,00000 

0.00000 

o.nonon 

U.OOOOO 

0,00000 

o.ooono 

..j. 06385 

0 . 057117 

0.25177 

11.37221 

i-,349e>* 

0.20234 

0.00000 

O.OOOOO 

o.oooon 

o.oooon 

o.oonoo 

0.00000 

.0,16239 

-0.05446 

0,16097 

0.34929 

V.3541 7 

u. 21129 

0.00000 

0.00000 

O.oooon 

o.ooOoo 

O.Onooo 

0.40000 

-0.11785 

-0, 05661 

0.ii920(> 

u. 20234 

0,21129 

•>.12736 


